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CHAPTER  1 


0 


I 


INTRODUCTION 

In  the  early  70’ s  minimax  criteria  for  frequency  domain 
design  of  digital  filters  received  a  vast  amount  of  attention 
([1],  [2]).  The  most  popular  among  these  formulations  of  the 
problem  is  the  one  due  to  Parks  and  McClellan  [3],  probably 
because  it  allows  a  very  large  class  of  design  specifications 
and  it  is  available  in  an  efficient  implementation  [4] . 

parks  and  McClellan  characterized  the  problem  of  the 
design  of  linear-phase  finite  impulse  response  (FIR)  digital 
filters  as  the  following  Chebyshev  minimization  problem: 

rain  max  ^W(eJ  )|  D(eJ  )  -  H(eJ  )1  i  (1.1) 

th(i)Jia()  <d4} 

where 

W(eJ  )  is  a  weight  function  (real  and  positive), 


D(e 


J  )  is  the  desired  objective  function, 


v-1 

-  1  -jiLl 

hfl)e  ,  where  M  is  the  order  of  ch#» 

1*0 


filter, 


3i  - 


i-1 


where  3. 

l  3 


<  *  <  i. 


[0,H; 


SI 


SI 


2 


and  Bj  •• 


*  j  *  i 


B.  j  ■  i 

i 


i,j  ■  1  )2 )  •  •  <L 


with  L  denoting  the  number  of  passbands  and  stopbands. 


A  i  ou 

Fig.  1  gives  a  pictorial  illustration  of  the  criterion.  W(e  ) 

.  »  i  U)  -  id)  . 

weights  the  error  |  D(e  )  -  H(e  ) |  on  5.  Parks  and  McClellan 

a  i  CD 

mainly  conceive  piece-wise  constant  behavior  for  W(e  ),  although 
they  leave  open  the  possibility  of  other  (strictly  positive) 

A 

behaviors.  Ihe  reason  why  W(eJ  )  cannot  take  on  the  value  0 

A  j(JU 

will  be  explained  in  Chapter  2,  D(e  )  is  a  piece-wise  constant 
function  for  the  case  of  multiple  passband-s topband  filters. 

A 

3  is  taken  so  that  it  doesn’t  contain  any  discontinuity  points 

A  jU> 

of  D(e  ).  The  phase  linearity  of  the  frequency  response  of  the 

a  4  *1) 

filter  H(e  ),  implies  chat 


H(eJ  )  =  Q(eJ  )  H(eJ  )  eJ 


(1.2) 


j  j  (1) 

where  3  is  constant  and  both  Q(e  )  and  H(e  )  are  real  functions 

f  ' 

defined  [i]  in  the  following  way :  v  7 


'In  (1.2),  for  -  real  constant,  the  phase  term  is  actually  e-  ' "  i- 
we  neglect  to  take  into  account  phase  jumps  of  size  2*  as  iz  is 
customarv  in  the  literature. 


( 


Case  2:  N  even 


itu 

Q(eJ  )  -  cos  '■» 


M-l 


H(eJ  )  =»  2  a(K)  cos(Kx) 

K=0 


where  M  ■  -|-and  a(K)'s  are  defined  by 


h(M-l)  -  -J-  a (0)  +  i  a(l) 


h(M-K)  =»  -5-  a(K-l)  +  -j-  a(K) 


K*2,3, . . . ,M-1 


h(0)  =  y  a (M-l) 


Substituting  (1.2)  into  (1.1)  results  into  the  expression  given 
in  (1.3) 

r  »  V<*»  It 

Iniax  W(eJ  )|0(aJ  )  -  H'eJ  )ji 


min 


(1.3) 


M 


la'K)J  n  'U'-i 
K.-U 


wnere 


\  X  nUJ 


W(e  )  =  W(eJ  )  Q(e  ) 


i'J) 


0(.J“)  -  ^ 

Q(eJ  ) 


5  ij  i<0|w«  J  and  Q(eJ  )  0i 


The  sets  ,  i  =  1,2,...,L  defined  above  are  called  pass 


i  X  1  X 

bands  if  D(eJ  )  =  i,  ±$3 1  or  3topbands  if  D(e-'  )  =  0,  '^€3^.  The 


secs  !]  x  <  x  <x  j,  with  i  =*  1,2,...,  L-l  constitute 

1  '  si  ft 


the  so  called  "don't  care"  bands. 
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Equations  (1.1)  and  (1.3)  have  the  structure  of  a  minimax 
(or  Chebyshev)  approximation  problem  for  generalised  (trigonometric) 
polynomials.  The  solution  to  this  class  of  problems  is  character¬ 
ized  by  the  Alternation  Theorem  ([5]  ,  pp.  75'  E.  Ya  Remez  [6] 
provided  algorithms  for  the  computational  S(  tion  of  such  Chebyshev 
problems.  J.  McClellan  et  al.  published  a  c  uter  program  [4] 
that  uses  the  2nd  Remez  Exchange  Algorithm  ..nd  the  solution 

of  a  discretization  of  (1.3).  Specifically,  the  problem  solved 
in  McClellan's  program  is: 

min  max  bl(e^  )|D(e^  )  -  H(e^  )|J  (1.4) 

.  M  uje<i 

ia(K)i.=0  d 

with  £  [a-jujeJ  and  ill  3  u  m  -  -  ,  K  positive  integer) 

C  is  a  parameter  (positive  integers)  that  can  be  controlled  by 
the  user.  It  can  be  shown  ([5],  pp.  89-95)  that  the  solution  to 
(1.4)  approaches  the  solution  to  (1.3)  as  a  tends  to  0.  McClellan 
suggests  that  values  of  L  —  16  produce  satisfactory  results. 

The  choice  of  solving  (1.4)  instead  of  (1.3)  is  dictated  for 
reasons  of  computational  efficiency. 

McClellan's  program  works  very  well  for  designing  low-pass 
or  high-pass  digital  filters.  However,  the  multiband  filters, 
i.e.,  the  filters  with  a  number  of  pass  and  stop-bands  greater 
than  2,  designed  with  McClellan's  program  often  exhibit  non 
mono  conic  '  sometimes  even  resonant")  behavior  in  the  "don't 
care"  bands.  Figures  7>  11,  15,  19  illustrate  several  cases 


of  this  phenomenon. 


Since  McClellan's  program  became  the  tool  most  universally 
used  for  the  design  of  finite  impulse  response  filters,  Rabiner, 
Shafer  and  Kaiser  [7]  addressed  the  multiband  design  problem  in 
particular.  The  technique  they  propose  is  the  following:  if 
McClellan's  program  returns  a  multiband  design  with  resonances, 
che  filter  specifications  passed  to  the  program  are  modified 
according  to  empirical  strategies  until  a  filter  without  resonances 
is  obtained.  Their  strategies  take  into  consideration  the  modi¬ 
fication  of  the  size  of  the  stopbands  (i.e.,  basically  changes 
of  <f)  as  well  as  the  modification  of  W(e  ).  The  number  of  tries, 
with  the  McClellan's  program,  r^ecessary  to  obtain  an  acceptable 
filter  varies  from  case  to  case.  The  implementation  of 
multiband  filters  with  this  procedure  might  easily  become 
cumbersome . 

The  present  work  reconsiders  the  problem  of  the  minitnax 
design  in  the  frequency  domain  of  Linear  phase  finite  impulse 
response  (FIR)  digital  filters.  The  objective  is  to  provide  a 
satisfactory  theoretical  solution  for  the  design  of  multiband 
filters  as  well  as  a  convenient  technique  for  the  implementation 
v?  f  such  a  sc  Lu cion. 

It  was  decided  to  use  McClellan’s  program  for  the  imple¬ 
mentation  of  the  new  soLution.  This  was  natural  since  the  new 
algorithm  is  an  extension  of  the  Parks  McClellan  technique  and 


because  the  new  algorithm  would  be  of  interest  to  the  great  number 
of  filter  designers  currently  using  McClellan's  program.  Chapter  2 
shows  that  the  inadequacy  of  McClellan's  program  for  the  design  of 
multiband  filters  lies  in  the  formulation  of  the  problem.  Chapter  3 
presents  a  new  formulation  capable  of  handling  these  difficulties. 
Chapter  4  introduces  an  implementation  of  the  new  solution  and 
examines  the  results.  Chapter  5  points  out  the  relationship 
between  the  new  program  and  the  program  CONRIP  [8] . 
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CHAPTER  2 

LIMITS  OF  PARKS  AND  McCLELLAN'S  FORMULATION 


Two  standard  results  of  Approximation  Theory  are  now 

nd 

introduced:  the  Alternation  Theorem  and  the  2  Remez  Algorithm. 

They  respectively  constitute  the  theoretical  and  the  computational 
tools  to  solve  Chebyshev  problems. 

Alternation  Theorem:  Let  5  be  any  closed  subset 
or  [0,tt]  and  let  D(e^)  be  any  continuous  real  valued  function 
defined  on  3.  In  order  that  H*(eJ  )  be  the  unique 
best  approximant  on  3  to  D(e^X),  among  the  class  of  the 

trigonometric  polynomials  of  order  M,  it  is  necessary  and  sufficient 

iaj 

that  E(e  ),  defined  as 

E (e^X)  =  W(ej,JJ)|D(eJ'U)  -  H*(ejtU)|  (2.1) 


exhibits  on  3  at  least  M  +  2  "a  ltemations" .  Thus 


rju 


L\  x-lN  .  -  •;  .max  |_ ,  ja\  I 

H(e  )  38  -t(e  )  *  +  E  *:  =  +  |E(eJ  )| 

with  Xq  £  £  •••  5  and 

Proof:  The  proof  can  be  found  in  [5],  pp.  75. 

Remark:  The  notation  of  the  following  chapters  is 

consistent  with  the  one  of  the  previous  sections.  3,  H,  D,  W 
therefore  are  as  defined  in  (1.3). 

Notice  that  the  Alternation  Theorem  is  an  existance 
theorem,  it  does  not  describe  how  to  find  the  best  aoproximant 


ft 


9 


1  Cl) 

H*(e  ).  Remez  algorithms  are  iterative  procedures  for  generating 

i  t-\ 

H*(e  ).  The  second  one  applies  to  classes  of  approximating 
functions  like  the  trigonometric  polynomials. 

2nd  Remez  Algorithm:  Each  step  of  the  algorithm  works  with  a 
set  of  M  +  2  frequencies  The  frequencies  are  arbitrarily 

chosen  at  the  first  step  and  updated  at  successive  iterations 
according  to  the  particular  algorithm  ([5],  pp.  97).  The  frequencies 
^K^K=*0  are  use<*  to  determine  the  following  system  of  M  +  2  equations 
in  the  M  +  1  a(K)  's  and  in  p  (M  +  2  unknowns): 


E(ejli)K)=w(ejx'K) 


DCe-^)  -  H(e^K) 


(-1)  P 


(2.2) 


for  K  =  0,1,2,...,  M  +  1 


The  assumption  W(eJ  )  >  0,  ulcJ  allows  (2. 2)to  be  written  as: 

d 


1  COSOj  cos2ojq  ....  COSMO)q 


cos 


"1 


cos2. 


i  •• 


.  .  cosMu;. 


.•1+1 


,>i+1....cosMuM+1 


W(eJX0) 


W(eJxl) 


(-1) 


M+l 


V(eJ  'jf+1) 


a(0) 


a(l) 


I  a(M) 
i  i 


D(ejx0) 


i Wl 

D(eJ  ) 


’  W“M) 


D(eJ “M+l) i 


j  - 


(2.3) 


10 


For  reasons  of  efficiency  Parks  and  McClellan  avoid  the 
matrix  inversion  in  the  solution  of  (2.3).  This  is  done  by  first 


calculating: 


I**1 

2  bj  D(eJ  ) 

1»Q  _ 

(-D1  bt 

i=0  WCe-J^i) 


(2.4) 


K  ^  1 

where  b  =•  (-1)  /  - — — — 

K  i=0  xi  *K 


=»  COS,iJ^ 


Then  the  Lagrange  interpolation  formula  in  the  baricentric 

j  U)  r  -  M 

fora  is  used  to  interpolate  H(e  )  on  the  M  +  l  points  lx J 

K  KSU 

to  obtain  the  values 


CR  -  D(ejajK)  -  (-1)K 


W(ej% 


(2.5) 


with  K  *  0,1,... ,M. 

This  type  of  interpolation  gives  an  equiripple  fit  to  the  M  + 


data  points.  This  fora  of  interpolation  is  required  by  the  2nd 


Remez  algorithm. 
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K 


E 


■ 

E 


After  the  system  is  solved,  the  algorithm  calls  for  the 

check; 

E(ejJ,d)  <  E(ejU^)  -  p  ,  <*»d«3d  (2.6) 

If  there  is  some  frequency  of  3^  not  satisfying  inequality  (2.6), 
a  further  iteration  that  begins  with  the  update  of  the  frequencies 

is  necessary.  The  algorithm  continues  until  inequality 

-  i  **» 

(2.6)  is  satisfied  over  the  whole  <Jj.  The  H(e  )  obtained  at 

jU) 

this  step  is  the  best  approximant  H*(e  )  characterized  by  the 
Alternation  Theorem. 

The  Parks  and  McClellan's  FIR  design  technique  corresponds 
to  the  application  of  the  Alternation  Theorem  and  of  the  2nd  Remez 
algorithm  to  the  solution  of  (1.3). 

It  must  be  noticed  that  the  "don't  care"  bands  constitute 
a  mathematically  ambiguous  feature  of  their  formulation  that  can 
be  very  dangerous.  In  fact,  one  actually  does  "care"  about  the 

iUJ 

behavior  of  H(eJ  )  over  the  whole  band  [0,  ■■]  and,  in  particular, 

VJU 

one  wants  H(e  )  to  be  monotonic  over  the  "don’t  care11  bands. 

Further,  it  is  not  clear  that  ignoring  the  behavior  of  the 
filter  over  the  "don’t  care"  bands  will  necessarily  result  in  the 
desired  mono tonic  response. 

Before  showing  in  detail  the  dangers  connected  with  the 
"don’t  care"  bands,  let’s  review  a  few  mathematical  concepts 
necessary  to  understand  them.  Recall  that  the  extrema  of  a  functio 
over  a  compact  and  bounded  set  can  only  occur  at  the  boundary  point 
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of  the  set  or  at  the  interior  points  where  the  first  derivative 

i  u* 

of  the  function  is  0.  H(eJ  )  is  a  trigonometric  polynomial  of 
order  M  defined  on  [O,17].  Its  derivative  ^  can  have  at 

most  M  -  1  distinct  0's  in  (O,77)  since  it  is  a  polynomial  of 
order  M  -  1.  The  boundary  points  of  M(e^  )  are  the  set  lo,77). 
Therefore  H(eJ  )  can  have  at  most  M  +  1  distinct  extrema  on  [0,77] . 
The  function: 


Jw- 


.  i  cu  V 

E(eJ  )  =  H(eJ  )  -  D(eJ  ) 


(2.7) 


i  *  ~ 

where  D(eJ  )  and  are  defined  in  formulation  (1.3), is  a  trigonometric 

j  0) 

polynomial  of  order  at  most  M  on  J.  Notice  that  E(e  )  over  [C,  tt] 
is  not  a  polynomial  of  order  M,  because  of  the  discontinuities  of 

iU) 

D(eJ  ). 

The  cases  of  the  low  and  high-pass  filters  will  be  con¬ 
sidered  first  because  of  their  special  characteristics.  Then  the 
multiband  filters  will  be  discussed. 

For  low  or  high-pass  filters  «?  »  [0,oj  ]  ‘J  [  w  ,"] ,  the 

o  s 

interior  of  is  ^  h  (0,0^  )^(  ws>:t)  a^d  boundary  of  is 

2,  .  ;  \ 

o'  s  10,^  E(e  )  can  have  at  most  M  -  1  distinct  extrema 

®  "  j  JU 

on  «J*  (since  E(eJ  )  is  a  polynomial  of  degree  at  most  M  in  <*)  and 
at  most  four  extrema  on  the  boundary  1 0,  *',Jg .  Therefore 


E(e  )  can  have  at  most  M  +  3  extrema  on  Furthermore  (2,7) 

j  x 

implies  that  all  the  extremal  points  of  Z(e  ),  except  possibly 

j  * 

the  extrema  at  and  ^  are  extremal  points  of  H(e  ),  namely 
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-  H.  (.J1")  -  D(«J“) 


«3  (2.8) 


will  have  at  least  M  +  2  distinct  extrema  alternating  in  3  by 
the  Alternation  Theorem.  Specifically,  E*(e  )  can  have  at 
least  M  +  2  distinct  extrema  on  3  either  having  M  -  1  distinct 

O 

extrema  in  3  and  at  least  3  extrema  in  ,“>,*}  or  having 

0  P 

M  -  2  distinct  extrema  in  3  and  four  extrema  inlo,<" 

P  3 

4(1) 

More  specifically  E*(e  )  can  have: 


i)  M-l  extrema  in  3  and  extrema  at  Co,*11  ,if'j 

P  s 

ii)  M-l  extrema  in  3  and  extrema  at 

iii)  M-l  extrema  in  3  and  extrema  at  to,'xs>ir} 

iv)  M-l  extrema  in  3  and  extrema  at 

v)  M-l  extrema  in  3  and  extrema  at  ^p'^s 

vi)  M-2  extrema  in  3  and  extrema  at  lo^p'^g'*}. 

Some  of  the  cases  of  E*(eJ  )  listed  above  leave  open  the  possi 
bility  of  a  corresponding  H*(e  )  that  would  be  unacceptable 
as  a  low  or  high-pass  filter.  Case  ii) ,  for  instance  could 

■jU) 

correspond  to  and  H*(e  )  having  an  extremum  at  *  (see  Fig.  2 

Case  iv)  could  bring  an  H*(e^  )  with  a  saddle  point  in  (<*>  , <*>  ) 

P  s 

JOB 

(see  Fig.  2f)  and  case  vi)  an  H*(e  )  with  a  local  maximum  in 
(^pj'-^g)  (see  Fig.  2e) ,  xhe  reasons  such  possibilities  don't 
occur  are  explained  in  the  following  theorem. 


t 
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Theorem:  The  low  and  high-pass  digital  filters  designed  via 

the  Remez  algorithm,  according  to  formulation  (1.3)  where 

B  £  <  u;  <  a)  ,«w  € £0 , t ] }  is  the  don't  care  band,  have 

P  s  P  s 

the  following  properties: 

a)  and  ug  are  always  extremal  points  of  E*(eJ  ), 

(i.e.,  the  cases  ii)  and  iii)  above  cannot  occur) 

i  uj 

b)  and  oj  are  not  extremal  points  of  H*(eJ  ) 

jui 

c)  H*(e  is  strictly  monotonic  on  B,  namely: 
dH*(eju:) 

- — -  <  0  ,  U)€B  for  the  low-pass  filters 

Q,  Ui 

iOU 

dH*(eJ  ) 

and  - dui -  **  0  >  ^B  for  the  high-pass  filters 

Proof:  the  Proof  takes  into  consideration  each  of  the 

juu 

cases  of  E*(e  )  listed  above  and  shows  the  necessity  of 
properties  a),  b),  c)  for  the  cases  whose  occurrence  is 
not  a  contradiction  (i.e.,  all  but  case  ii)  and  iii)). 


Case  i)  satisfies  properties  a),  b~>,  c) 

Property  a)  is  satisfied  by  definition. 

i  Hi  ^  i  X 

E*(eJ  )  has  now  M  +  3  extrema  on  J.  Therefore  H*(eJ  ) 


has  >1  +  1  distinct  extrema  in 
and  2  at  0  and  ”,  respectively, 
therefore  cannot  be  extrema  of  H* 


namely  M  -  1  in  5 


The  points  x  and  * 
P 

J  ^ 

(a  )  (property  b)) 


and 


further  there  cannot  be  any  extremal  point  in  15  (property  c)). 
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Case  11)  and  111)  contradict  the  Alternation  Theorem 

Consider  first  case  ii) ,  i.e.  EMe^)  has  M-l  extrema 

O 

on  &  plus  3  extreme  at  {0,  u>  , it  } .  Notice  that  0  and  it  are 

P 

also  extrema  of  H*(e^u),  therefore  — L  -  0  cannot  occur 

d  a) 

at  any  tu  e  B  (property  c)).  Also,  the  point  m  cannot  be  an 

P 

extremal  point  of  H*(e^u)  otherwise  H*(e^u)  also  has  M  +  2 

extrema  on  5  .  Since  ,  by  assumption  <u  is  not  an  extremal 

s 

point  of  E*(e^U)»  E*^111)  does  not  alternate  (see  fig.  2a, b) 
and  therefore  the  contradiction  is  achieved. 

A  completely  analogous  argument  proves  that  case  ill)  cannot 


occur. 


Case  iv)  and  v)  satisfy  properties  a) ,  b) ,  c) 

1 00 

Consider  case  iv)  first,  i.e.  E*(e  )  has  M-l  extrema  in 
plus  3  extrema  at  '<•0,  <^,  u  }.  Notice  that  0  is  extremal  point 
also  of  H* (e^ J).  therefore  H*(e^)  has  at  least  M  extrema  in 


It  is  worthy  to  distinguish  3  subcases: 


Subcase  I:  -  and  uj  are  both  also  extrema  of  H*(e^").  This  is  a 
-  p  s 

contradiction  because  H*(e^)  would  have  M+2  extrema. 

Subcase  II:  is  an  extremum  of  H*(e^)  and  r  is  not  (the  case 

-  P  s 

w  is  an  extremum  of  H* (e^ and  a  is  not, is  completely  analo- 
s  P 

gous  and  leads  Co  the  same  conclusions  )• 

*  .  0  ^ 

3v  hypothesis  H*(eJ  "0  has  M-l  extrema  on  "  ,  since  H*(eJ  ) 

trerr.a  on  (0,  ~)  it  is  dHA(e^  “)/d-  ^  0  W  ^  :  5 


has  at  most  M-l  ex 


\' 


(recall  that  the  roots  of  derivatives  correspond  to  the  interior 
extremal  points). 


*!• 


W8)*H  O  (mf3)*3 


FIs;.  2a, b.  Case  ii):  with  M  -i  extrema  in  ~ 

and  three  extrema  at  ro,xr>,r}.  The  frequency 
o3  is  not  extremal  point* of  E*(e-^). 


In  this  subcase  the  contradiction  is  achieved  because 


E*(e^)  violates  the  Alternation  Theorem  (see  Fig.  2-c,d). 
Therefore  subcase  II  cannot  occur. 

Subcase  III:  Neither  to  nor  u  are  extrema  of  H* (e^) .  In 
-  p  s 

this  case  property  a)  and  b)  are  satisfied  by  assumption. 

H*(e^w)  by  assumption  has  M-l  extrema  in  3  >  therefore  jff* } 

do) 

can  not  occur  in  B  otherwise  H*(e^w)  would  have  M  extrema 
in  (0,  t) .  This  proves  property  c) . 

The  same  argument  shows  also  that  the  properties  a) ,  b) , 
c)  are  satisfied  for  case  v)  (it  is  enough  to  interchange  the 
roles  of  0  and  ir  as  extremal  and  non-extremal  points  of  EMe^) 
respectively) . 


Case  vi)  satisfies  properties  a),  b) ,  c) 


E*(e^w)  is  assumed  to  have  four  extrema  at  I  0, 


%  l  '  p-  S' 

and  M-2  extrema  in  Ir  .  This  implies  that  H* (e^u)  has  2  extrema 


at  0  and  tt  plus  M-2  extrema  in?,  i.e.  M  extrema  in 

It  is  again  convenient  to  distinguish  three  subcases. 


T  u: 

Subcase  I:  ^  and  ^  are  extrema  of  H*(e  ).  This  cannot  occur 
s  p 


because  H*(e^  )  would  have  M+2  extrema  in?. 


Subcase  II:  ^  is  an  extremum  of  H*(e^  )  and  is  not  (eauiv- 
p  s 


alent  to  the  case  ^  is  an  extremum  of  H*(eJ  )  and  -  is  not) 

S  D 


H*(eJ~)  has  by  assumption  M  -  2  extrema  in  ~ ,  two  extrema  at 


’*0,  "  .v  and  one  extremum  at  -  .  H*(e~  )  therefore  has  bv  assumn- 

? 


tion  M  ^  1  extrema  in 


This  excludes  the  possibility  of 


E*^)  o  H*(e'") 


0 

Fig.  2c, d.  Case  iv):  E^Ce^)  with  M  -1  extrema  in  c 

and  three  extrema  at  {o,u)p,o)s}.  The  frequency 
yjp  is  assumed  to  be  extremal  point  of 
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—  — L  _  n  at  some  u  e  B  (property  c) ,  otherwise  H*(e^W) 

du» 

would  have  M  +  2  extrema  in  3  .  In  this  subcase  contradiction 

is  achieved  because  E*(e^U)  violates  the  Alternation  Theorem 

(the  situation  is  analogous  to  the  one  of  Fig.  2-c,d). 

Subcase  III:  u>  and  u  are  not  extrema  of  H*(e^U).  This  case 
-  p  s 

specializes  into  four  situations: 

(a)  Maximum  (or  minimum)  in  B.  I.  e.: 


dH(ejtl>) 


d2H(ejai) 


#  0  ,  oj  e  B 


(2.8) 


Since  H*(eJ  P)  j  H*(eJ  s)  as  long  as  oj  ^  oj  and  neither  oj  nor 

p  s  p 

u>s  are  extrema  of  'H*(e^u>),  the  presence  of  a  maximum  (or  a 

minimum)  in(uj  ,  oj  )  implies  also  the  presence  of  a  minimum 
p  s 

(or  a  maximum)  in  (u  ,  u>  )  (see  Fig.  2e).  The  contradiction 

P  s 

is  achieved  since  this  requires  H*  (e~^ u)  to  have  M  +  2  extrema 


on  [  0  ,tt  ]  , 

(b)  Saddle  point  in  B  (see  Fig.  2f) 


dH*(ej,Jl)  d2H*(ejai) 

du  .2 

d  oj 


(J  £  3 


(2.9) 


Recall  that  a  saddle  point  of  a  polynomial  H(eJ  )  corresponds 

to  a  zero  of  the  derivative  of  multiplicity  at  least  2.  By 

i  u)  * 

assumption  H*(e  )  has  M  -  2  extrema  in  this  means  that 

1  jL)  ® 

dH*(e"  )/du  has  M  -  2  zeros  in  J  .  Condition  (2.9)  implies  that 
dH*(eJ  )/d,j  has  a  zero  at  least  of  multiplicity  2  in  B.  This 
is  a  contradiction  because  it  implies  that  dH*(e^u‘) /d(j  has  M 


zeros  in  (0,~)  and  dH*(eJ^)/duj  is  a  polynomial  of  order  M  -  1. 


H*(en 


0 

Fig.  2e.  Case  vi) :  E*(eJw)  with  M  -2  extrema  in  U 

and  four  extrema  at  {o,oip,u)s,}  H*^01)*3 

assumed  to  have  a  maximum  in  ((*)„, u)c) . 

P  ® 


Fig.  2f.  Case  vi) :  E'Ce^)  with  M  -2  extrema  in  c 
and  four  extrema  at  r o ,:js ,t}  H*(e^)  is 
assumed  to  have  a  saddle-point  in  (^p,us). 


9 


(c)  Further  extremum  of  H* (e^ W)  in  7,  i.e.  if  c 

O 

denote  the  set  of  the  extremal  frequencies  of  3  assumed 
hypothesis,  this  circumstance  is  expressed  as 


by 


dH*(ejui) 

du 


0  at  u  e 


to  e  *s 


M-2 

u  i  •  }  . 

i-1 


This  situation  satisfies  properties  a),  b) ,  c)  by  assumption. 

It  should  be  noticed  that,  since  w  is  not  accounted  among 

M“  2 

the  extremal  frequencies  { cu  ,  the  Alternation  Theorem 

1  i-l 


imposes 


E*(eJui)|  <  p 


max  |E*(e^u) 


|E(ejui) 


,  M-2 


(d)  No  further  extremum  on  (0,  it  ).  I.e.,  if  {a  }*  denotes 

o  1  i”l 

the  set  of  the  extremal  frequencies  of  7 ,  this  circumstance  is 
expressed  as 


M-2 

- - ;  i  0  £  |  0J  £  ( 0,^0  ,  f  ^  ,  WjKWg,  Up} 


This  situation  satisfies  properties  a) ,  b) ,  c)  by  assumption. 


Remark:  For  simplicity  E(e^w)  has  been  used  in  the  theorem  according 

to  definition  (2.8)  instead  of  definition  (2.1).  The  theorem  can  be 

S  (jj 

proved  also  for  a  weighted  error  as  in  (2.1)  with  W(e  )>  0  and  real 

Also  in  this  case  Ete^)  turns  out  to  have  all  its  extremal  points, 

except  possibly  m  and  jj  ,  in  common  with  H(e^~)  and  E*(e^)  has 
s  p 

extrema  as  in  the  six  cases  previously  considered. 


H*(en 


Q 


B 


Fig.  2g.  Case  vi) :  with  M  -2  extrema  in 

and  four  extrema  at  {o,ajp,o3s,ir}  H* (e^)  ;is 
assumed  to  have  a  further  extremum  in  Sr  . 


t 


» 


o 
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The  preceding  theorem  shows  that  the  mathematical  structure  of 
the  approximation  problem  (1.3)  corresponding  to  the  design  of  high 
and  low-pass  filters  guarantees  strictly  monotonic  behavior  in  the 
"don't  care"  band.  It  is  legitimate  to  ask  if  formulation  (1.3) 
also  guarantees  strict  monotonicity  in  the  "don't  care"  bands  for 
multiband  filters.  The  answer  is  "no".  The  reasons  behind  it  can 
be  illustrated  by  an  example.  Consider  a  3-band  filter  like  the  one 


shown  in  Figure  3-a,b.  In  this  case  3 


-[o.  i 


the  interior  of  3  is  <5  *  (0>  ^  ^**f  us2^  ^  ^wf3*  an<^ 
boundary  of  7  *  (0,  Ms2'  UJf3’  11 

E*(^J:J)  given  by  the  2nd  Remez  algorithm  has  at  least  M  +  2 
extrema  in  7.  Therefore  E*  can  have  the  structure  of  any  of  the 
cases  below: 


(1) 

M 

- 

4 

extrema 

in 

3 

and 

6 

extrema 

in 

^5: 

1 

subcase 

(2) 

M 

- 

3 

extrema 

in 

o 

3 

and 

5 

extrema 

in 

$3  : 

6 

subcases 

(3) 

M 

- 

3 

extrema 

in 

o 

3 

and 

6 

extrema 

in 

1 

subcase 

o 

(4) 

M 

- 

0 

extrema 

in 

and 

4 

extrema 

in 

'  J  : 

13 

subcases 

(5) 

M 

- 

o 

extrema 

in 

0 

.t 

and 

5 

extrema 

in 

0 

c  J  . 

6 

subcases 

o 

.  0 

(6) 

M 

- 

2 

extrema 

in 

and 

6 

extrema 

in 

tV  rs 

*j  . 

1 

subcase 

(7) 

M 

- 

i 

extrema 

in 

0 

and 

3 

extrema 

in 

0 

'J'  tT  . 

20 

1  subcases 

0 

9 

(3) 

M 

- 

i 

extrema 

in 

c 

and 

4 

extrema 

in 

r* 

15 

i  subcases 

(9) 

M 

- 

i 

extrema 

in 

0 

'  J 

and 

5 

extrema 

in 

O 

r. 

6 

subcases 

(10) 

M 

— 

i 

extrema 

in 

0 

and 

6 

extrema 

in 

O 

,-V 

1 

subcase . 

» 
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Each  one  of  these  cases  can  specialize  in  many  subcases  depending 
from  the  points  of  S'  7  that  are  assumed  to  be  extrema  (similarly  to 
what  was  done  for  the  low  or  high-pass  filters).  The  number  of  the 
subcases  is  6!/(K!  (6-K)  !) ,  where  K  is  the  number  of  extrema  in  3^3. 

It  is  not  difficult  in  this  forest  of  cases  to  find  a  situation 
unable  to  guarantee  monotonic  behavior  in  the' don't  care' bands  without 
violating  the  Alternation  Theorem  or  the  relationships  between  the 
order  of  a  polynomial  and  the  number  of  its  extrema.  Select,  for  exam— 
pie,  the  subcase  of  case  9)  corresponding  to  M-l  extrema  in  2  and 

extrema  at  (0,  u)g^,  ufl^’  usl'  uf2'  us2’  ^3  are  not  extrema 

of  HMe^10)  has  M  extrema  in  7  :  one  at  0  and  M-l  in  J. 

H*(eJU))  can  have  a  further  extremum  in  *  as  shown  in  Fig*  3  a,b 

without  causing  any  contradiction. 

In  general  for  an  L-band  filter 


cl 


L 

U 

j-1 


r 


and  E(e^)  can  have  up  to  M-l  +  2L  distinct  extrema  in  ~  .  The  2nd 
Remez  algorithm  guarantees  only  that  E*(eJuJ)  has  at  least  M  +  2 
alternating  extrema  in  5 .  The  possibility  that  the  Remez  algorithm 
takes  into  account  extrema  of  EMe-^)  not  corresponding  to  extremal 
point  of  H*(e^J)  is  very  high.  Therefore  for  an  L-band  filter  (L  >  2) 
some  of  the  M  +  1  extrema  of  H(e''")  can  occur  anywhere  on[0,’  ].  If 
they  occur  in  ~  the  second  Remez  algorithm  constrains  them  to  give 
deviations  from  D(eJ'L)  within  the  final  minimax  error  p  =  max  E*.(e'*'). 


If  they  occur  in  the"don't  care"bands  they  are  unconstrained  and  can 
not  only  spoil  the  monotonicity  but  also  become  resonances  for  the 
filter.  Unbounded  extrema  of  H*(e^U)  in  the"don't  care "bands  are 
reported  [7]  to  occur  experimentally  in  9  out  of  10  multiband  filters 
designed  with  formulation  given  in  (1.3). 

In  the  next  section  a  new  formulation  of  the  filter  design 
problem  is  presented.  This  new  formulation  is  immune  from  the  draw¬ 
backs  of  the  McClellan  formulation  and  capable  of  the  straightforward 
design  of  the  multiband  filters. 


CHAPTER  3 


NEW  FORMULATION  OF  THE  LINEAR-PHASE  FIR  DIGITAL  FILTERS  DESIGN  PROBLEM 


The  elimination  of  the  "don't  care"  bands  in  the  formulation  of 
the  linear  phase  FIR  digital  filter  design  problem  calls  for  the 
following  type  of  formulation: 


min  max  WCe^)  |  D ( e^ )  -  H(e^“)  | 

{a(k)}^Q  we  [0,  tt  1 


(3.1) 


where  the  symbols  are  as  defined  in  (1.3). 

The  presence  of  D(e^^),  as  defined  for  (1.3),  i.e.,  as  a  piece- 
wise  constant  discontinuous  function,  brings  two  major  difficulties 
to  the  formulation  (3.1).  The  first  difficulty  is  that  the  discon¬ 
tinuities  of  D^^)  will  seriously  limit  the  convergence  of  the 
approximation  (3.1).  For  instance  it  is  easy  to  see  that  if  W(e^W)*l 
the  minimax  error  will  not  become  smaller  than 


max 

jj  £  discontinuity 


TlD(ejU)I. 


(3.2) 


*1'j0  \ 

points  of  D(e  )> 


The  second  difficulty  is  that  the  Alternation  Theorem  cannot  be 
used  to  characterize  the  solution  H*^10)  of  (3.1),  since  it  requires 
the  continuity  of  the  function  to  be  approximated.  It  will  be  noticed 
that  since  the  approximating  functions  in  (3.1)  are  trigonometric 
polynomials  it  does  not  even  make  sense  to  require  a  discontinuous 
behavior  from  them.  These  reasons  motivate  the  use  of  a  continuous 
u(e^^)  in  (3.1).  Figure  ^a  shows  an  example  of  a  piecewise 


4a  -  Example  of  DCe^)  for  the  McClellan's  formulation 

b  -  Example  of  continuous  Dfe^)  for  the  new  formulation 
(technique  L) 

c  -  Example  of  continuous  DCeJ^)  for  the  new  formulation 
(technique  P) 
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D 

ft 


c 
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discontinuous  D(e^U)  and  Fig.  4b  and  c  two  continuous  versions  of 
DCe^10)  for  use  in  the  new  formulation  (3.1).  A  method  for  choosing 
such  a  continuous  DCe^)  is  another  major  point  of  the  new  formula¬ 
tion  of  the  filter  design  problem,  whose  detailed  discussion  will  be 
given  later  in  the  section. 

The  relationship  between  the  Parks  and  McClellan's  formulation 
(1.3)  and  the  new  formulation  will  now  be  discussed.  Such  a  compar¬ 
ison  turns  out  to  be  rather  informative  and  gives  the  motivation  for 
the  choice  of  W(e^"°)  and  D(e^u)  in  the  new  formulation  (3.1) 

The  Parks  and  McClellan  formulation  (1.3)  can  be  obtained  as  a 
particular  case  of  the  new  formulation  (3.1).  This  equivalence  occurs 
when 

4 

W(eJai)  -  0  u  eD  -  {u|ooe[0,  *3,  u*  3  }  (3.3) 

is  used  in  (3.1). 

It  will  now  be  shown  why  the  formulation  of  (1.3)  is  mathemat¬ 
ically  preferable  to  the  formulation  (3.1)  with  condition  (3.3), 
although  the  two  are  absolutely  equivalent  in  meaning. 

Formulation  (3.1)  for  W(e^U)  >  0  constitutes  a  minimization 
problem  with  respect  to  a  weighted  minimax  norm  defined  for  real 
functions  supported  by  [0,  t“(*)  Condition  (3.3)  turns  (3.1)  into 

(*)  A  norm  is  a  functional  £  over  a  vector  space  X,  with  the  following 
properties : 

(i)  f(x)  >  0  for  all  x  e  X 

(ii)  f(ax)  =  ;  a :  x  for  all  scalars  i  and  :c  s  X 

( iii)  f (x-*-y)  _  f (sc)  +  f(y)  for  each  x.y  t  X 

(iv)  £(x)  =  0  if  and  only  if  =  0. 
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a  minimization  problem  with  respect  to  a  weighted  seminorm  defined 
for  the  real  functions  supported  by[o;ir  ] .  Such  a  seminorm  defines 
a  norm  for  the  Co,  ir ]  real  functions  supDorted  by  3(**).  To  write 
(3.1)  with  condition  (3.3)  as  formulation  (1.3)  corresponds  to  con¬ 
sidering  the  seminorm  for  function  of  [0,  ir  ]  as  a  norm  for  functions 
of  T5.  This  point  of  view  is  of  theoretical  as  well  as  computational 
utility.  The  theoretical  utility  comes  from  the  fact  that  formula¬ 
tion  (1.3)  entitles  us  to  apply  to  a  seminorm  problem  the  results  given 
for  the  norm  problems  (such  as  the  Alternation  Theorem  and  the  second 
Remez  algorithm) .  In  order  to  understand  the  computational  utility 
notice  that  W(e^w)  *  0  for  u>  e  and  considered  as  a"don't  care" 
band  will  produce  the  same  effect,  namely  eliminate  the  occurrence  of 
extremal  frequencies  on  B^.  In  fact,  W(e^w)  m  0  will  force  E(e^w)  m  0 
on  j  e  B^,  therefore  the  second  Remez  Algorithm  will  not  find  any 
extremum  of  E(eJ  )  for  ut  c  B^.  If  B^  is  a  "don't  care"  band,  it  is  just 
not  taken  into  account  during  the  operation  of  the  second  Remez  Algo¬ 
rithm,  therefore  it  cannot  deliver  any  extrema  of  E(e^u).  Therefore 
the  total  effect  of  the  two  techniques  is  identical,  but  the  use  of 
the  "don't  care"  band  is  more  efficient.  The  efficiency  lies  in  saving 
the  actual  evaluation  of  E(e^)  over  u  £  B^  as  well  as  in  allowing 
the  efficient  non-conventional  solution  of  system  (2.3)  via  (2.4)  and 
(2.3)  (made  possible  by  the  fact  that  the  weight  W(e^'“)  is  never  0  in  3). 

(**)  A  seminorm  is  a  functional  g  over  a  vector  space  X,  satisfying 
the  properties  (i),  (ii),  (iii)  above. 


A  computational  check  of  the  above  claimed  equivalence  between 
"don't  care"  bands  and  regions  supporting  0  weight  was  tried. 

McClellan's  program  was  used  because  minor  modifications  can 
convert  it  into  a  tool  implementing  formulation  (3.1)  with  condition 
(3.3).  A  direct  verification  of  the  equivalence  was  found  not  to  be 
possible  since  the  McClellan's  program  does  not  allow  0-weights.  This 
limitation  derives  from  the  calculation  of  p  via  (2.4).  Also  an  indirect 
verification  by  means  of  small  weights,  ideally  tending  to  0,  was  not 
easy  to  obtain.  Specifically,  filters  of  relatively  high  order  (above 
60)  can  call  for  weights  of  order  10  ^  or  less  in  order  not  to  exhibit 
any  extremum  over  the  "don't  care"  bands.  McClelland  program  starts 
to  lose  its  numerical  accuracy  for  weights  of  this  order (as  reported 
in  [7]  )  and  the  results  may  not  be  reliable.  However  for  filters  of 
lower  order  the  indirect  verification  is  possible,  as  the  examples 
of  Fig  5  and  6  show. 

During  these  experiments  it  was  noticed  that  a  complete  removal 
of  the  extremal  frequencies  from  the  transition  regions  gives  a  better 
performance  in  terms  of  reducing  ripple  over  the  "care"  bands  than  a 
partial  removal  of  extremal  frequencies.  It  was  noticed,  in  practice, 
that  a  fewer  number  of  extremal  frequencies  in  the  transition  regions 
results  in  a  smaller  ripple. 

The  idea  behind  the  choice  of  the  values  of  W(e^~)  and  D(e^) 
over  the  transition  regions  for  use  in  (3.1)  is  to  have  as  few  extremal 
frequencies  as  possible  over  these  regions  (possibly  none).  The  choice 
of  the  values  of  r.v(e^)  and  D(eJ  )  over  the  pass-bands  and  the  stop- 
bands  follows  the  usual  criteria  taken  for  the  Parks  and  McClelland 
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formulation. 

The  W(e^)  to  be  assigned  to  the  transition  regions  comes  from 
a  trade-off  between  two  conflicting  requirements.  W(ejaJ)  should  be 
small  in  order  to  keep  E(e^U)  small,  so  that  no  extremal  frequency 
is  detected  by  the  Remez  algorithm.  But  WCe^10)  cannot  be  too  small, 
since  the  smaller  W(e^w)  is,  the  larger  the  term  iDCe^)  -  H(e^u)  | 
becomes  for  a  given  minimax  error  E(e^).  As  a  matter  of  fact,  the 
difficulty  with  W(e^)  =  0  (or  "don't  care"  bands)  is  that  resonances 
of  HCe^)  (i.e.  points  where  the  term  |  D (e^  -  H ( e^ ^ )  |  is  very  large) 

can  occur  without  being  detected  as  extrema  of  E(e^).  Therefore  the 
weight  over  the  transition  regions  has  to  be  taken  small,  but  non-zero, 
in  order  to  prevent  resonances. 

The  following  two  observations  are  useful  for  the  choice  of 
D(e^)  over  the  transition  regions.  The  first  observation  is  that 
every  multiband  filter  can  be  thought  to  be  the  composition  of  several 
low-pass  or  high-pass  filters  [  7  }  Such  low  and  high-pass  filters  will 
be  called  "prototypes"  in  the  present  work.  The  second  observation  is 
that  the  prototypes,  by  the  theorem  of  Chapter  2  have  a  monotonic 
transition  region.  Thus,  they  can  be  safely  implemented  with  the 
McClellan's  procran. This  suggests  that  the  transition  regions  of  the 
prototypes  be  taken  as  a  model  for  the  transition  regions  of  D(e^). 

Two  techniques  for  obtaining  the  transition  regions  of  D(e^~)  have  been 
tried*  a  muitisegment  piece-wise  linear  approximation  of  the  transition 
regions  of  the  prototvpes,  and  the  direct  use  of  the  transition  regions 
of  the  prototypes  in  Die""").  For  convenience  these  will  be  rererred 
to  as  technique  L  and  ?,  respectively .  A  detailed  discussion  or  their 


implementation  and  performance  will  be  given  in  the  next  section. 

The  new  formulation  together,  with  technique  P  prevents  an  extre¬ 
mal  frequency  from  occurring  in  the  transition  region  for  the  design 
of  low  or  high-pass  filters,  since  the  error  E(e^)  over  the  trans¬ 
ition  region  is  made  0  by  the  term  jDCe^)  -  HCe^)).  This  result  is 
independent  of  the  value  of  W(e^)  over  the  transition  region. 

In  the  design  of  multiband  filters,  the  complete  elimination  of 
the  extremal  frequencies  from  the  transition  regions  is  not  to  be 
expected,  even  with  technique  P,  since  every  prototype  induces  extremal 
frequencies  into  the  transition  regions  of  the  other  prototypes. 

Finally,  notice  that  the  new  formulation  (3.1)  allows  a  general 
control  over  the  transition  regions  of  the  filter.  This  characteristic 
can  be  used  to  obtain  monotonic  behavior  as  well  as  any  other  desired 
behavior  in  the  transition  regions.  The  new  formulation  is  therefore 
very  suited  to  the  design  of  filters  for  which  the  shape  of  the  trans¬ 
ition  regions  is  important.  Thus  if  transition  band  performance  is 
important  in  a  high  or  low-pass  filter  formulation,  (3.1)  is  to  be 
preferred  to  the  McClelland  formulation.  Figure  7  shows  an  example 
of  a  low-pass  filter  that  could  not  be  obtained  with  McClelland  form- 
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CHAPTER  4 

IMPLEMENTATION  AND  PERFORMANCE  OF  THE  NEW  FORMULATION 

The  new  formulation  of  the  linear  phase  digital  filter  design 
problem  (3.1)  has  the  structure  of  a  minimization  problem  in  minimax 
norm  for  trigonometric  polynomials,  similar  to  the  Parks  and  McClellan 
formulation  (1.3).  The  Alternation  Theorem  and  the  second  Remez 
algorithm  are  still  applicable  for  computing  the  solution  to  problem 
(3.1). 

The  program  written  by  J.  McClellan  contains  a  general  purpose 
subroutine  for  performing  the  Remez  algorithm.  The  program  is  designed 
in  three  parts;  an  input  section,  a  computational  part,  and  an  output 
section.  The  first  part  is  devoted  to  building  W(e'510),  D(e^w)  from  the 
input  data  and  to  setting  up  the  approximation  problem.  The  central 
computational  part  is  the  implementation  of  the  Remez  algorithm.  The 
third  part  is  devoted  to  the  display  of  the  results.  Such  modularity 
facilitates  possible  modifications  of  the  program. 

It  appeared  convenient  to  incorporate  the  implementation  of  the 
new  formulation  (3.1)  into  McClellan's  program.  This  would  make  the 
new  formulation  directly  accessible  to  *.1  who  currently  use  McClellan's 
program. 

The  key  idea  of  the  implementat ion  was  to  insert  in  parallel 
with  McClellan's  program  an  alternate  pattern  in  3  parts  devoted  to 
the  implementation  of  (3.1).  The  first  part  can  accept  the  input  data 
peculiar  to  the  new  formulation.  The  second  part  prompts  the  operation 
of  the  second  Remez  Algorithm  over  the  full  band  [0,  -  The  third 
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part  is  devoted  to  the  display  of  the  results  of  interest.  The  new 
program  so  obtained  can  be  thought  of  as  a  modified  version  of  the 
McClellan's  program  allowing  the  selection  of  two  different  "modes" 
of  operation.  Namely  the  "regular  mode"  that  prompts  the  new  program 
to  implement  formulation  (1.3)  (i.e.  to  behave  exactly  as  the 
McClellan's  program) and  the  "custom  mode"  that  prompts  the  implementa¬ 
tion  of  the  new  formulation  (3.1).  A  user-oriented  description  of  this 
program  is  presented  in  Appendix  1. 

The  rest  of  the  section  is  devoted  to  the  presentation  of  the 
results  obtained  with  the  new  program  together  with  some  practical 
observations  useful  for  the  choice  of  W(e^u)  and  D^^)  over  the  trans¬ 
ition  regions. 

For  comparison  purposes  the  four  filters  reported  in  [7]  as 
typical  cases  of  multiband  filters  with  non-raonotonic  transition 
regions  have  been  designed  with  the  McClellan '  s  program, the  new  program, 
and  CONRIP  (which  is  .mother  digital  filter  design  program  discussed 
in  Chapter  5).  The  filters  examined  are  labeled  Design  1,  2,  3,  or 
4  as  in  [7  }  The  new  program  has  been  used  with  both  the  techniques 
L  and  P,  described  in  Chapter  3. 

Figures  8-23  constitute  by  themselves  the  best  comments  on  the 
performance  of  the  new  program  versus  the  other  two  programs.  It 
can  be  seen  that  the  new  program  gives  strictly  monotonic  behavior  in 
the  transition  regions  also  in  cases  where  McClellan's  program  and 
CONRI?  do  not. 

Several  comments  are  now  presented  about  the  choice  of  W<’eJ~) 
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Fig.  10a  New  program,  technique  P: 

H*(e^U>)  of  Design  1 
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Fig.  10b  New  program,  cechnique  P:  extremal 

frequencies  of  Design  1 


Fig.  10c  New  program,  technique  P: 

H*^^)  of  Design  1  in  dB 
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Fig.  11 


CONRIP:  H*(ejx)  of  Design  1 
(from  [8] ,  pp.  66) 


3.2 
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3.5 
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ig.  12b  McClellan's  program:  extremal 

frequencies  of  Design  2 
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Fig.  18a  New  program,  technique  P: 

iU) 

H*(eJ  )  of  Design  3 
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Fig.  13b  New  program,  technique  ?:  e>:tremaL 

frequencies  of  Design  3 


Fig.  20b 


McClellan's  program:  extremaL 
frequencies  of  Design  4 
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Fig. 


!lc  New  program,  technique  P: 

H* (e^  )  of  Design  4  in  dB 
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New  program,  technique  P:  desired 

i  ^ 

function  D(eJ  )  for  Design  4 
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and  DCe^)  for  the  new  program  when  technique  L,  i.e.,  the  multisegment 
piecewise  linear  simulation  of  the  transition  regions,  is  used. 
Technique  L  implies  the  subdivision  of  each  transition  region  into 
several  bands  over  which  D(e^W)  has  a  prescribed  slope.  The  values 
of  DCe"^)  on  each  band  are  chosen  to  linearly  approximate  the  filter 
prototypes.  A  criterion  for  the  choice  of  the  weights  on  each  band, 
that  seems  very  effective,  is  to  assign  the  smallest  weights  to  the 
internal  sub-bands  of  the  transition  region  where  is  steepest. 

Greater  weights  should  be  assigned  to  the  other  sub-regions  where  the 
absolute  value  of  the  slope  of  D(e^U)  become  smaller.  For  instance, 
if  the  transition  region  between  the  band  with  weight  ot  and  the  band 
with  weight  2  is  divided  into  5  sub-bands  the  preceding  criterion 
can  be  expressed  by  the  two  following  choices  of  the  weights  (it  is 
assumed  :t  >  3,  without  any  loss  of  generality). 

Choice  1  Choice  2 

Band  1  a  ':- 


Subregion  1 

set 

0.  t 

Subregion  2 

ssa 

O 

Q 

Subregion  3 

sssa 

0.00a 

Subregion  4 

ssS 

O.Oi 

Subregion  5 

si 

O.i 

3and  2 

Choice  2  is  a  particular  case  of  Choice  1  and  it  has  been  extensively 
used  in  the  design  examples  presented  in  this  section. 
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Choice  1  assumes  a  >  sa  >  ssa  >  sssa  and  ss3  <  sS  <  S  and  it  is 
recommended  when  Choice  1  gives  poor  results  (for  instance  in  Deisgn 
4).  A  non-uniform  (piecewise)  constant  choice  of  the  weights  in  the 
transition  regions  seems  to  significantly  contribute  to  the  elimina¬ 
tion  of  extremal  frequencies  on  them.  Technique  P  corresponds  to 
assigning  the  values  of  the  transition  regions  of  the  prototypes  to 
the  corresponding  regions  of  D(e^).  The  weight  assigned  to  a  trans¬ 
ition  region  with  this  technique  should  be  the  smallest  weight  still 
capable  of  giving  monotonic  behavior. 

The  order  of  the  prototypes  is  a  point  that  needs  some  comment. 

If  the  multiband  is  thought  of  as  a  cascade  of  prototypes,  their  order 
should  be  equal  to  the  order  of  the  multiband  divided  by  the  number 
of  the  prototypes.  If  instead  the  multiband  is  thought  of  as  a  parallel 
of  prototypes  their  order  should  be  equal  to  the  order  of  the  multi¬ 
band.  Actually  the  relationship  between  the  multiband  and  the  proto¬ 
types  is  not  clear.  Therefore,  the  question  of  the  order  of  the  proto¬ 
types  doesn't  have  a  definite  answer.  Fortunately,  it  appears  that 
the  order  of  the  prototypes  doesn't  influence  the  result  very  much, 
probably  because  of  the  effect  of  the  weight  that  helps  keep  E(eJ^) 
snail  over  the  transition  regions.  In  the  design  examples  shown  in 
this  section  the  prototypes  were  usually  taken  of  the  same  order  of 
the  multiband.  In  Design  4,  though,  the  multiband  is  of  order  73 
and  the  prototypes  were  taken  of  order  41.  Prototypes  of  different 
orders,  according  to  the  empirical  rules  of  Rabiner  et  al.[  7  ]  were 
tried,  but  no  improvement  over  the  other  choices  cf  the  order  was 


not  iced . 
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The  plots  of  this  section  show  that  technique  P  gives  better 
monotonicity  of  the  transition  regions  than  technique  L.  However, 
the  two  techniques  have  comparable  performance  in  terms  of  elimina¬ 
tion  of  extremal  frequencies  from  the  transition  regions.  This  is 
explained  by  the  two  following  facts:  every  prototype  induces  extremal 
frequencies  into  the  transition  regions  of  the  other  prototypes, 
therefore,  as  mentioned  earlier,  extremal  frequencies  in  the  transi¬ 
tion  regions  should  be  expected.  Furthermore,  technique  L  benefits 
from  the  non-uniform  weights  in  the  transition  regions,  while  tech¬ 
nique  P,  in  the  present  implementation,  doesnTt  have  this  feature. 

(This  further  modification  has  not  been  implemented  because  the  results 
of  technique  P  were  already  satisfactory). 

To  obtain  a  multiband  filter  with  monotonic  transition  regions 
with  the  new  program  is  very  simple.  Some  common  sense  is  needed  to 
choose  the  weights  of  the  transition  regions,  which  is  the  only  euristic 
part  of  the  procedure.  The  rules  are  the  following:  if  the  initial 
weights  give  resonances,  then  increase  their  value;  if  the  initial 
weights  give  monotonic  behavior  then  try  a  smaller  weight  that  might 
reduce  the  ripple.  The  criteria  for  the  choice  of  the  weights  in  the 
transition  regions  should  be  used  in  these  changes  of  weight. 

All  the  filter  examples  shown  in  this  section  were  obtained  on 
the  first  try,  with  the  exception  of  the  filter  of  Design  i  which 
required  two  attempts  (it  is  not  excluded  that  prototypes  of  order  71, 
instead  of  ^1,  would  have  riven  a  sat isf actorv  result  on  the  first  trv). 
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CHAPTER  5 

COMMENTS  ABOUT  CONRIP 


This  section  presents  the  relationships  between  the  new  program 
discussed  in  Section  4  and  the  program  CONRIP,  written  by  M.  T. 
McCallig  [8],  CONRIP  implements  the  design  of  FIR  filters  according 
to  a  "constrained  ripple"  formulation,  due  to  B.  J.  Leon  and  M.  T. 
McCallig  [8].  Their  formulation  is  the  following:  given  continuous 
functions  U(e^w)  and  L(e^U)  on  [  0,  it  ]  such  that  U(e^)  >  LCe^)  find 
the  polynomial 

.u)  M 

P(e^  )  =  Z  a(K)  cos(Kcj)  such  that: 

K=0 

(i)  L(ej“)  <;  P(ejuJ)  <;  U(ejw),  u>  e  [0,  tt  ] 

(ii)  P(eJUJ)  is  monotone  on  specified  subintervals  of  [0,  tt] 

(iii)  P(e^J)  is  the  minimal  order  polynomial  meeting  conditions  (i) 
and  (ii). 


It  is  worthwhile  to  note  that  the  preceding  formulation  is  not  a 
traditional  approximation  problem  (like  (1.3)  and  (3.1))  and  it  is  a 
full  band  formulation  (like  (3.1)  and  unlike  (1.3)). 

CONRIP  allows-  the  design  of  multiband  filters  that  do  not  exhibi 
strictly  mcnotonic  behavior  in  the  transition  regions  but  which  dc 
not  have  resonances  like  the  ones  obtained  with  McClellan’s  program. 
The  results  of  using  CONRIP  for  the  filters  of  Designs  1,  2,  3,  1  are 
shewn  in  figures  11,  15,  19,  and  23. 


'he  similar  it" 


ripple  problem  and  the 


between  the  constrained 
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Chebychev  approximation  problems  (1.3)  and  (3.1)  is  greater  than  it 
might  appear  at  first  sight.  In  fact  the  computational  solution  of 
the  constrained  ripple  problem  is  obtained  with  an  algorithm  (reminis¬ 
cent  of  the  second  Remez  algorithm)  which  searches  for  the  polynomials 
tangent  to  U(e^u)  or  Lte^)  at  their  extrema.  Two  theorems  of 
McCallig  state  precisely  the  similarity  of  the  two  methods. 

Theorem  5.1  Let  H*(e^u)  be  the  solution  to  the  Chebychev  approx¬ 
imation  problem  (3.1).  Let 


p  =  max 

U)£  Co.tt] 


W(eJj  I  H* (eJ  )  -  D(eJ  ) 


(5.1) 


The  constrained  ripple  problem  having 


L(eJai)  =  D(ej“)  -  p 
U(e^w)  =  D(ejuJ)  +  p 


(5.2) 


admits  unique  solution  P(e'^UJ)  =  H*(e^U). 

Proof:  [3],  pp.  42.  The  converse  of  Theorem  5.1  stated  as 
follows . 

Theorem  5.2:  Let  the  boundary  curves  U(e^J)  and  L(e^)  be  given 
such  that  the  unique  solution  to  the  constrained  ripple  problem 
is  ?(ejAJ).  The  Chebychev  approximation  problem  (3.1),  having 


D(ej")  =  4(U(ej'J)  +  L(ej*)) 


W(eJ ")  = 


U(ej°)  -  D(eJ°) 
U(eJ~)  -  Q(eJ~) 


admits  unicue  solution  K*(eJ~)  = 


i« 


Note:  The  assumption  U(e^  )  >  L(e^  )  guarantees  the  denominator 

U(ej  )  -  D(ej  )  =■  h  (U(ej  )  -  L(ej  ))  >  0. 

Proof:  [8],  pp.  42. 

It  is  important  to  notice  that  Theorem  5.1  says  that  the  approx¬ 
imation  problem  (3.1)  needs  to  be  solved  in  order  to  obtain  the  equiv¬ 
alent  constrained  ripple  problem.  (The  two  problems  are  said  to  be 
equivalent  if  they  have  the  same  solution).  In  fact  the  final  devi¬ 
ation  p  (5.1)  is  not  available  before  the  solution  of  (3.1). 

Theorem  5.2  similarly  says  that  the  constrained  ripple  problem 
needs  to  be  solved  in  order  to  obtain  the  equivalent  approximation 
problem.  The  order  of  the  polynomial  to  use  in  (3.1)  is  not  avail¬ 
able  without  solving  the  constrained  ripple  problem.  Some  consequences 
of  the  above  two  theorems  having  practical  interest  are  now  c mussed. 

The  following  equalities 


L(ejw) 

U(eJw) 


D(ejui) 

D(ej*) 


W(e^w) 


W(e^) 


(5.4) 


can  be  used  to  solve  via  CONRIP  the  following  Chebvchev  approximation 
problem : 


min  min  max  T.v(eJ~)  1  D(e^’°)  -  H(e^)  1  (5.5) 

M:z"  •.«) 

This  last  aspect  of  the  equivalence  between  constrained  ripple  prob¬ 


lems  and  approximation  problems  needs  a  comment. 
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Let 

p (0)  *  min  max  0W(e^w)  | Dte*^)  -  H(e^w)[  (5.6) 

{  ^k)}k-o  weC 0 ,Tr-* 

with  0  real  and  positive.  The  best  approximant  H*(eJW)  of  (5.6) 
is  independent  of  0,  but  the  minimax  error  p(0)  is  dependent  on  0. 

The  constrained  ripple  formulation,  as  intuitively  understood  and  as 
precisely  stated  in  Theorem  5.1,  is  sensitive  to  the  minimax  error 
p(S).  Therefore  the  filters  found  via  formulation  (3.1)  and  the  new 
program  are  independent  of  the  weights.  However,  scaling  the  weights 
in  formulation  (5.5)  implemented  via  CONRIP  by  (5.4),  give  different 
filters  (even  differences  in  the  order  of  the  filter  must  be  expected). 

The  new  formulation  (3.1),  on  the  converse,  suggests  different 
ways  of  using  CONRIP.  Perfect  monotonicity  of  the  transition  regions 
should  be  obtained  by  picking  L(e~^)  and  iKe^10)  via  (5.3)  where  the 
DCe^)  is  chosen  with  the  help  of  the  prototypes  and  the  W(eJ^)  is 
chosen  with  the  criterion  of  Section  4.  The  new  program  presented 
in  Section  4  can  be  used  to  implement  the  following  problem  which  is 
related  to  a  constrained  ripple  problem: 

Given  continuous  functions  U(e^~)  and  L(e^)  on  [  0,~]  such  that 

U(e^)  >  L(e^),  and  a  given  order  M,  find  the  polynomial 
M 

?(e^)  «  I  a(K)  cos  (fw)  such  that: 

K=0 

(1)  L(e^)  -:(M)  -  P(e:")  <  U(e:~)  +  :(M)  for  some 

SOI)  >0  ^  e  [  0,  -  ] 

(ii)  ?(e"^)  is  monotone  on  specified  subintervals  of  [0,"]. 


It  should  be  noted  that  the  6(M)  of  point  (i)  depends  on  M, 
and  there  exists  an  M  such  that  y  M  >  M  6(M)  *  0. 

The  experimental  verification  of  the  use  of  CONRIP  to  solve 
approximation  problem  (5.5)  and  of  the  use  of  the  new  program  to  solve 
problems  similar  to  the  constrained  ripple  problem  was  beyond  the 
objective  of  this  work.  Nevertheless,  such  verification  would  be 


an  interesting  project. 
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CHAPTER  6 
SUMMARY 

The  problem  addressed  in  this  work  was  finding  a  technique  for 
the  straightforward  design  in  the  frequency  domain,  using  a  minimax 
criterion  of  multiband  linear  phase  FIR  digital  filters. 

To  the  author’s  knowledge  the  only  similar  attempt  has  been  the 
work  of  Rabiner,  Kaiser  and  Schaffer  [7].  Their  work  has  the  merit 
of  having  brought  to  general  attention  the  problems  arising  in  multi¬ 
band  FIR  filter  design  with  McClellan’s  program.  However,  their 
solution  doesn’t  address  the  essence  of  the  problem  which  lies  in  the 
inadequacy  of  Parks  and  McClellan’s  theoretical  formulation  of  the 
filter  design  problem  for  the  multiband  filter  case.  The  ’’strategies” 
proposed  in  [7  ] to  design  multiband  filters  are  empirical  and  their 
application  is  generally  not  straightforward. 

Chapter  2  presents  an  original  analysis  of  Parks  and  McClellan’s 
filter  design  formulation.  It  is  clearly  stated  that  its  mathematical 
meaning  over  the  band  [0,:r  1  is  that  of  an  approximation  problem  not 
according  to  the  Chebvchev  norm  but  according  to  a  seminorm.  The  Parks 
and  McClellan  formulation  is  shown  to  be  ideal  from  the  filter  design 
point  of  view  for  the  high  and  low-pass  filter  cases.  It  is  similarly 
shown  that  it  is  mathematically  inadequate  from  the  filter  design 
point  of  view  where  there  is  more  than  one  ’’don’t  care”  band  as  in  the 
multiband  filter  case. 

Chapter  3  proposes  a  new  formulation  of  the  filter  design 
problem.  The  new  formulation  corresponds  to  a  minimization  problem 


» 


L 


in  Chebychev  norm  over  the  full  band  [0,  ].  It  requires  the  use  of 

continuous  desired  functions  DCe^).  The  Alternation  Theorem  and  the 
second  Remez  algorithm  still  apply.  The  new  formulation  is  immune 
from  the  drawbacks  of  the  McClellan  formulation  and  it  is  theoretically 
adequate  for  the  multiband  filter  design  problem. 

Chapter  4  discusses  an  implementation  of  the  new  formulation 
based  on  McClellan’s  program.  Considerations  about  the  choice  of 
W(e'JUJ)  and  D(e^)  having  practical  interest  are  also  introduced.  The  per¬ 
formance  of  the  new  program  is  compared  to  that  of  McClellan's  program  and 
McCallig's  program  CONRIP  for  the  multiband  filter  case.  It  is  empha¬ 
sized  that  the  new  program  seems  to  be  the  only  one  capable  of  giving 
strictly  monotonic  transition  regions  in  a  straightforward  way  without 
changing  filter  specifications. 

Chapter  5  presents  the  relationships  with  McCallig's  program 
CONRIP.  The  possibility  of  the  use  of  CONRIP  to  implement  a  design 
criterion  very  close  to  the  one  presented  in  Chapter  3  and  the  possi¬ 
bility  of  the  use  of  the  new  program  to  implement  a  design  criterion 
similar  to  a  constrained  ripple  problem  are  both  introduced  and  dis¬ 
cussed. 

From  the  author’s  point  of  view  McClellan's  filter  design  formu¬ 
lation  (1.3)  is  a  particular  case  of  the  new  formulation  (3.1).  This 
made  it  natural  to  incorporate  McClellan's  program  as  the  core  of  a 
program  implementing  formulation  (3.1).  The  practical  result  obtained 
in  this  thesis  therefore  is  a  kind  of  "extended  McClellan's  program" 
that  is  identical  to  the  original  one  when  appropriate,  such  as  for 


the  low  and  high-pass  filters,  or  can  be  used  for  its  new  features 
(formulation  (3.1))  when  McClellan's  program  is  not  adequate,  as  in 
the  design  of  multiband  filters. 


70 


i  , 

APPENDIX  1 

a 

Appendix  1  provides  a  user-oriented  description  of  the  new 
program  discussed  in  Section  4.  The  program  operates  interactively 
I!  and  asks  the  user  a  sequence  of  questions  part  of  which  are  derived 

from  McClellan's  program  and  part  are  original.  The  meaning  of  the 
questions  in  terms  of  McClellan's  program  will  not  be  reviewed  here 
tm  and  the  word  ‘'standard”  will  be  used  in  place  of  the  actual  answer 

for  questions  from  McClellan's  program.  All  the  comments  to  the  answers 
will  be  on  the  modifications  for  the  new  program.  The  questions  will 
be  numbered  for  convenience  of  reference. 

1.  TYPE  FILTER  ORDER 
Standard 

n 

2.  ENTER  FILTER  TYPE 
Standard 

3.  ENTER  NUMBER  OF  BANDS 

Enter  the  number  of  bands  where  DCe^)  changes  slope  if  the  transi¬ 
tion  regions  will  be  piecewise  linearly  simulated.  Enter  the  total 
number  of  passband,  stopband,  and  transition  bands  if  the  transition 
regions  of  the  prototypes  will  be  used. 

4.  TYPE  GRID  DENSITY 
Standard 

3.  ENTER  IPRINT  IPLCT 
Standard 

t  6.  ENTER  LINE  PRINT  FLAG  (0  -  DO  NOT  PRINT,  1  -  PRINT) 

Standard 


* 
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7.  STANDARD  WEIGHT  AND  TARGET,  TYPE  0,  CUSTOM  WEIGHT  AND  TARGET,  TYPE  1 
0  selects  the  operation  of  the  program  as  regular  program  of 
McClellan,  1  selects  the  "custom  mode"  operation,  implementing 
formulation  (3.1)  (all  the  comments  to  the  questions  assume  1  is 
entered  here) 

8.  TO  CHANGE  BAND  EDGES  TYPE  1,  OTHERWISE  TYPE  0 

The  band  edges  previously  entered  will  be  kept  if  0  is  entered 
(this  happens  because  the  program  can  cyclically  call  itself). 

New  band  edges  must  be  provided  if  1  is  entered. 

9.  ENTER  THE  BAND  EDGES  -  -  NO.  =  NBANDS 

Enter  the  sequence  of  the  edges  of  the  bands.  Since  the  new  pro¬ 
gram  assumes  full-band  operation  the  edges  corresponding  to  two 
contiguous  bands  have  to  be  1  sample  apart.  In  order  to  facilitate 
this  feature,  the  program  just  uses  the  edges  0,  0.5  and  the  even 
ones  (the  second  edge,  the  fourth,  and  so  on)  to  calculate  the  odd 
ones  via  an  increment  of  1  sample. 

Since  the  important  information  for  this  answer  is  the  edge  0, 
the  edge  0.5  and  the  even  edges,  the  odd  edges  are  usually  assigned 
to  dummy  integers  like  0  or  l(because  they  are  convenient  to  type!). 

10.  FILTERS  PROTOTYPE  USED?  NO  =  0,  YES  -  1 

Enter  1  for  technique  P.  Enter  0  if  no  filter  prototype  will  be 
used . 

11.  ENTER  SIMULATION  FLAGS 

Question  11  appears  only  if  1  has  been  entered  at  question  10. 

Enter  a  vector  associating  a  number  with  each  band.  If  the  band 
is  a  transition  region  the  number  associated  must  be  the  number  c: 
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the  file  containing  the  impulse  response  of  the  prototype  for  that 
region.  The  number  0  has  to  be  given  for  the  pass  and  stop-bands. 

12.  ENTER  ORDERS  OF  THE  PROTOTYPES 

This  question  appears  only  if  1  has  been  entered  at  question  10. 

Enter  a  vector  associating  a  number  with  each  band.  If  the  band  is 
a  transition  region  the  associated  number  must  be  the  order  of  the 
prototypes  used  for  it.  Any  number  can  be  associated  with  the  other 
bands.  This  feature  allows  the  use  of  prototypes  of  different  order 

13.  PIECEWISE  LINEAR  DESIRED  FUNCTION?  YES  *  1,  NO  *  0 

Enter  1  for  technique  L.  Enter  0  if  technique  L  is  not  wanted. 

It  should  be  pointed  out  that  the  current  implementation  allows 
also  the  use  of  technique  L  and  technique  P  together,  that  is  some 
transition  regions  can  be  taken  from  prototypes  and  some  others 
can  be  piecewise  linearly  modeled. 

14.  ENTER  DESIRED  FUNCTION  AT  THE  EDGES 

This  question  appears  only  if  1  has  been  entered  at  question  13. 

Enter  the  values  assumed  by  D(e^)  at  the  band  edges. 

15.  ENTER  THE  CONSTANT  VALUES  FOR  EACH  BAND 

This  question  appears  only  if  0  has  been  entered  at  question  13. 

Enter  1  corresponding  to  the  pass-bands,  0  corresponding  to  the  stop- 
bands  ,  and  any  number  corresponding  to  the  transition  regions, 
lb.  ENTER  WEIGHT  FACTORS  FOR  EACH  BAND 

Enter  the  weights  corresponding  to  each  band.  It  should  be  empha¬ 
sized  that  the  current  implementation  calls  for  uniform  weights 
ever  each  band.  A  piecewise-constant  weight  could  be  also  obtained 
over  a  region  by  means  of  its  subdivision  int?  several  bands. 
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Figure  24  shows  an  example  of  conversational  terminal  using  the  new 
program  with  technique  P. 
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Fig.  2 4-  Conversational  terminal  example 

using  the  new  program:  Design 
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APPENDIX  2 

This  appendix  offers  numerical  examples- of  filters  designed  with 
techniques  L  and  P.  The  examples  correspond  to  the  filters  denoted 
Design  1,  2,  3,  4. 

Both  techniques  require  one  to  design  filter  prototypes.  Tech¬ 
nique  L  needs  them  as  models  for  piecewise  linearly  simulating  the 
transition  regions.  Technique  P  needs  them  to  use  the  actual  values 
of  their  transition  regions  into  D(e^). 

The  prototypes  are  designed  with  McClellan’s  program.  Their 
specifications  are  directly  obtained  from  those  of  the  multiband  filter 
corresponding  to  them. 

For  every  design  example  the  information  relative  to  the  proto¬ 
types  wi i 1  be  presented  first,  and  then  the  information  relative  to 
the  multiband  filter. 

(a)  Technique  L 

Design  1 

Prototype  1.1  (low-pass);  Filter  order  =  75 


Band  no . 

First  Idee 

Second  Face 

D ~) 

W  ( e "  "  ) 

j 

0 . 165  2  394  7 

- 

■*  f  a 

0. 1 4  3  ~  3  ?  2 

0.5 
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Prototype  1.2  (high-pass)  ;  Filter  order  =  75 


1 

Band  no. 

First  Edge 

Second  Edge 

D(ej“) 

W(ej") 

1 

0 

0.41679744 

0 

0.03853275 

2 

0.37032451 

0.5 

1 

1 

The  transition  region  of  Prototype  1.1  was  modeled  by  means  of 

2  line  segments.  The  one  of  Prototype  1.2  was  modeled  by  means  of 

3  line  segments. 


Cl 
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Multiple  passband-stopband  filter  1 


Band  No. 

1st  Edge 

2nd  Edge 

D(ej'") 
(1st  Edge) 

D(ejw) 
(2nd  Edge) 

U’(eJ"  ) 

1 

0 

0.14375973 

1 

1 

0.04588809 

0 

0.159 

1 

0.07 

0.004593809 

3 

0 

0.165 

0.07 

0 

0 . 04533809 

4 

0 

0.37032451 

0 

0 

1 

5 

0 

0.  385 

0 

0.09 

0.038385327 

6 

0 

0.4021199 

0.09 

0.91 

0.003853275 

/ 

0 

0.41679744 

0.91 

1 

J- 

'  .03853275 

3 


0 


0.  3 


i 


0.0335:275 
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Design  3 

Prototype  3.1  (Low-pass);  Filter  order  *  57 

Band  No.  1st  Edge  2nd  Edge  D(e^W) 

1  0  0.11018835  1 

2  0.00820222  0.5  0 


Prototype  3.2  (High-pass);  Filter  order  =  57 

3and  No.  1st  Edge  2nd  Edge  D(e','°) 

1  0  0.26931373  0 

2  l).  2C585967  0.5  1 


E 


Prototype  3.3  (Low-pass);  Filter  order  =  57 

3and  No.  1st  Edge  2nd  Edge  ) 

1  0  0.37673551  1 

2  0. 31153715  0.5  0 


W(ejuJ) 

0.19386528 

0.17459027 


W(ej^) 

0.17459027 

1 


U*(eJ  ") 

1 

0.13180259 
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Prototype  3.4  (High-pass);  Filter  order  =  57 


3  and 

No. 

1st  Edge 

2nd  Edge 

D(ejul) 

W(ejU) 

1 

0 

0.46299174 

0 

0.18180259 

2 

0. 

3892995 

0.5 

1 

0.21319649 

The 

transition 

regions  of 

Prototypes  3.1, 

3.2,  3.3,  3.4  were 

simulated  by  3,3,5,  and  3  line  segments  respectively. 


Multiple  passband-stopband  filter  2 


Band  No, 

1st  Edge 

2nd  Edge 

D(ej'~) 
(1st  Edge) 

D(eJUi) 

(2nd  Edge) 

W(eJw) 

1 

0 

0.00820222 

1 

1 

0.19386528 

2 

0 

0.03439052 

1 

0.96 

0.0193 

3 

0 

0.04139057 

0.96 

0.91 

0.0193 

4 

0 

0.077 

0.91 

0.09 

0.00174 

5 

0 

0.084 

0.09 

0.04 

0.0174 

6 

0 

0. 11018835 

0.04 

0 

0.0174 

7 

0 

0.20585967 

0 

0 

0.17459027 

8 

0 

0.216 

0 

0.065 

0.017459027 

9 

0 

0. 2591734 

0.065 

0.935 

0.001749027 

10 

0 

0.26931373 

0.935 

1 

0.01749027 

2 1 

0 

0.31150715 

1 

1 

1 

i  •'> 

i.  _ 

0 

0. 31342266 

1 

0.953 

0.01 

13 

0 

0. 3699 

0.933 

0 .  047 

0.001 

14 

n 

0. 37673551 

0.047 

0 

0.018130233 

0 


0. 3391995 


0 


0. 13 130253 


td  No. 

1st  Edge 

2nd  Edge 

D(eJ  ) 
(1st  Edge) 

D(eJ  ) 

(2nd  Edge) 

W(eJ  ) 

16 

0 

0.4045 

0 

0.006 

0.018180258 

17 

0 

0.44779124 

0.006 

0.994 

0.0018180258 

18 

0 

0.46299174 

0.994 

1 

0.02139649 

19 

0 

0.5 

1 

1 

0.21319649 

(b)  Technique  P 
Design  2 

Prototype  2.1  (High-pass)  ;  Filter  order  =  27 


Band  No.  1st  Edge  2nd  Edge 

1  0  0.22754845 

2  0.0728033  0.5 


D(e^) 

0 

1 


W(ej“) 


0.1616032 


Prototype  2.2  (Low-pass)  ;  Filter  order  »  43 


Band  No.  1st  Edge  2nd  Edge 

1  0  0.3445328 

2  0.29030124  0.5 


D(eJUi) 

1 

0 


W(ejUJ) 
0. 1616032 
0.0034068 


The  prototypes  used  in  this  example  are  of  different  arders  (their 
orders  correspond  to  those  suggested  by  Rabiner  et  al.  [  7]). 

McClelland  program  stored  the  coefficients  of  Prototypes  2.1 
and  2.2  in  disk-files  number  21  and  22  respectively. 
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Multiple  passband-stopband  filter  2 

D(ej“) 

W(ejuJ) 

Band 

Simulation 

Prototype 

1st 

2nd 

No. 

Flags 

Orders 

Edge 

Edge 

1 

0 

0 

0 

0.07280333 

0 

1 

2 

21 

27 

0 

0.22754845 

1000 

0.01616032 

3 

0 

0 

0 

0.29030124 

1 

0.1616032 

4 

22 

43 

0 

0.3445328 

1000 

0.00340608 

5 

0 

■  0 

0 

0.5 

0 

0.00340608 

The  dummy  value  1000 

in  the 

transition  regions 

of  D(e^W 

)  has  been 

used  to  indicate  the  use  of  Prototypes. 


Design  4 


Prototype  4.1  (High-pass)  ;  Filter  order  -  41 


Band  No. 

1st  Edge 

2nd  Edge 

D(ejw) 

W(eju) 

1 

0 

0.13199438 

0 

0.0959953 

0 

0.08886197 

0.5 

1 

0.11187421 

Prototype  4.2  (Low- 

-pass)  ;  Filter 

order  =  41 

Band  No. 

1st  Edge 

2nd  Edge 

D(eJ~) 

W(eO 

1 

0 

0.27193968 

0 

0.11187421 

? 

0.18550831 

0.5 

1 

1 

Protoc} 

:pe  4.3  (High-pass)  ;  Filter 

order  =  41 

Band  No. 

1st  Edge 

2nd  Edge 

D(eJ ") 

W(eJ  “) 

1 

0 

0. 35373202 

0 

2 

0.28319105 

0.  5 

0. 11177379 
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Prototype  4.4  (Low-pass).;  Filter  order  ■  41 
Band  No.  1st  Edge  2nd  Edge  D(e^u) 

1  0  0.45732656  1 

2  0.43737502  0.5  0 


W(eja)) 

0.11177379 

0.05401694 


The  prototypes  in  this  case  were  taken  of  orders  different  from 
the  order  of  the  multiple  passband-stopband  filter.  McClellan’s 
program  wrote  the  coefficients  of  Prototypes  4.1,  4.2,  4.3,  4.4  in 
disk-files  number  21,  22,  23,  24  respectively. 


Multiple  passband-stopband  filter  4 


Band 

No. 

Simulation 

Flags 

Prototypes 
t  Order 

1st 

Edge 

2nd 

Edge 

D(eJ“) 

W(ejuJ) 

1 

0 

0 

0 

0.08886197 

0 

0.0959953 

o 

21 

41 

0 

0.13199438 

1000 

0.09 

3 

0 

0 

0 

0.18550831 

1 

0.11187421 

4 

22 

41 

0 

0.27193968 

1000 

0.11 

5 

0 

0 

0 

0.28819105 

0 

1 

6 

23 

41 

0 

0.35373202 

1000 

0.11 

7 

0 

0 

0 

0.43737502 

1 

0.11177379 

S 

24 

41 

0 

0.45732656 

1000 

0.05 

9 

0 

0 

0 

0.5 

0 

0. 05401695 

* 


APPENDIX  3 


This  Appendix  contains  the  computer  printouts  with  all  the 
numerical  information  referring  to  the  filters  shown  in  the  plots 
of  this  work-. 

In  order  to  facilitate  the  association  of  the  plots  with  the 
information  tables  corresponding  to  them,  the  tables  of  this  Appendix 
are  labeled  with  the  same  number  as  the  figures  containing  the  plot 
to  which  they  refer.  Therefore,  Figure  4  corresponds  to  Table  4, 

Fig.  13  corresponds  to  Table  13,  and  so  on. 

It  should  finally  be  noticed  that  the  printouts  of  the  new 
program  used  with  technique  L  under  the  voice  "DESIRED  VALUE"  report 
the  slopes  of  D(e^)  over  each  band. 
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TABLE  5 


'phtc  impulse  response  cptr) 
LI^Ea*  Pma3£  digital  ptlte*  design 
*eme i  exchange  algorithm 

BANQPAS3  PILtE« 

filter  length  •  3a 


•••#•  IMPULSE  RESPONSE 


*C 

n- 

• 

30) 

N( 

2)> 

a 

*? 

25) 

«C 

3)« 

a.3rjuj5«a«.oj 

• 

29) 

H( 

4)  • 

•9 •445o9436C«9a 

a 

N( 

27) 

Nf 

33* 

a 

"( 

2*) 

HC 

*)• 

a 

N( 

29) 

*c 

7)  ■ 

»,l3!*3!0»t.»l 

a 

Hf 

2«) 

H  c 

9)  ■ 

•9,l732Tt9«E.9t 

a 

N( 

23) 

N( 

*)• 

a 

22) 

*C 

1 0)« 

9.39»l7«l<e.91 

a 

*t 

at  3 

*c 

ill* 

a 

*< 

20) 

*C 

i2:« 

.9,3T33879«C.9t 

a 

19) 

*< 

133  • 

•9t«34T33l9e.<lt 

a 

"t 

19) 

*c 

U)« 

9,l4»77J<*fc««90 

a 

*■*( 

17) 

*c 

13)« 

9,«4948«37E*99 

a 

*< 

1*) 

d  a  no  l  dANO  2  9  anO 

L3*e»  3a.no  eqge  0,390000000  3.100009000 

UPPER  dANO  EDGE  3f2e?e3m9  0,500009903 

9ESIPE0  value  1,390030000  3,300000000 

WEIGHTING  1,300009900  1.303003003 

DEVIATION  a, 33436072a  3.3323%3T23 

OEVIaTIQN  in  oe  *53,71**205*5  *53.71**309*3 

E*7*r«AL  P*EQU£NCIES 


3,2000000 

0,3375000 

0,3750000 

0,1393331 

3,  I4t***7 

2,l7ad333 

a,  i*l***7 

a,  2300030 

3,3309000 

3.3393333 

3,32*16*7 

3,33*2509 

3,3975900 

0,4197300 

3,4590000 

3. 4633333 


TABLE  6 


•i*ir*  twitse  oeaFaNae  cfto) 

LI*€A*  F*aS£  OCSITal  FILt£9  58SI^ 
0g*£2  g*Cn4*G£  AL3C4 1 T hH 

SaN0^a33  FIL^C* 

FlL^C*  L£*Gt"  »  30 

•••••  i^uljc  ••••• 

*  C  i)«  -O.lJl5at04C.02  »  -*C  33) 

h<  ax  a.2i*»M2ie*42  ■  -(  29) 

«<  3) »  a#3«45«90d€*.aa  •  *c  20) 

*<  4)s  -0,a3l5j422£.02  •  Hf  27) 

h (  5)»  .a. 750930450.02  «  hc  2m 

H(  b)«  0,0900t429C.02  »  rtf  23) 

*c  7) ■  a, j4ab4a23£«dt  ■  *<  24) 

h(  ax  -0,x4353i0e-ot  •  m<  23) 

Hf  9) «  .a,2»49a2fet£-3l  •  H  (  22) 

*t  10)  •  0,294907090-01  a  «<  2!) 

*(  113*  0.43»l207a£.0i  •  H(  £3) 

k<  1 2) •  -a, 550599400.01  >  xj 
*t  i3)*  .*,a*7*323«e.ai  •  h(  i a ) 

h<  14)*  0,149247340*00  •  *C  t?) 

H(  13)*  3#4339€332£*33  •  "C  X) 

0ANO  l  04*0  2  34*0  3 

13*0*  34*0  0OG0  3.33333030V*  0,202003333  0,390003333 

JF»*S  9 4 NO  0OG0  0,230000800  3,330000009  3,300009030 

D03I90U  V  A  (.  J  £  0.300000909  -  t  3 , 300«U0000  3,339200000 

*£.G«r:»G  1.330000903  3,800010000  1,809909000 

C<TSr-4L  F0C^li£^Xg1 

3.3303909  3,8395333  ?,,)75C*00  0cl;04lft7  3,14X447 

3, 1 7a5  3  33  3,19X447  3.2300000  3,3020*33  8,3134X7 

0,3312590  *,  35*3333  3,3075900  a,M*75o0  8,4523633 

3,4833333 
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TABLE  7 


ri*jT£  t*PUL3€  9C3PQ n$E  CFt.9) 
L:ncaR  * n 4 Sg  o  .'JITaL  *TUTE»  3531G* 
iF.HE*  EXCHANGE  awgOPIThH 

9AN0PA53  PtLTEQ 

PlL'E®  LENGTH  m  40 


”(  1 3 •  *.i Jia;a3oe-J2 

iti 

a 

-C 

40) 

1( 

2)  * 

*,«49479S3E.03 

a 

*( 

39) 

*c 

3?  « 

-J.U413saU.3i 

a 

if 

30) 

1C 

•o » 

•9 , 5 1  6o9 1645-43 

a 

1C 

37) 

-c 

5 )  * 

-4.464550005-03 

a 

i  c 

3b) 

*c 

6  )  ■ 

•0,iT9b3a3lE-92 

a 

If 

35) 

*c 

7  5  « 

5,349250235.03 

a 

1C 

34) 

H( 

93* 

-0.15960056E-02 

a 

1C 

33) 

*c 

93  • 

0.1S7E2673E«02 

a 

1C 

32) 

HC 

133  • 

9.792S«290C.93 

a 

1  c 

31) 

*c 

U3  • 

0,11  6790325-92 

a 

H( 

30) 

1< 

1 23  • 

0.323t7824£.0a 

a 

1( 

29) 

*c 

1  33  * 

•9,31 4 l 3921E-92 

a 

1  ( 

20) 

1C 

143  * 

0,«48i^4t7C.a3 

a 

1( 

27) 

*c 

153- 

-0.T187A533E-92 

a 

HC 

26) 

1  ( 

lb)  ■ 

-0,135403625-91 

a 

1C 

25) 

*c 

173- 

0,296955235.42 

a 

If 

24) 

1( 

ld)  ■ 

-0,279652*55-01 

a 

1( 

23) 

*c 

191  ■ 

0,1  3a83039E*90 

a 

H( 

22) 

*C 

203  - 

0,433643235*30 

a 

1C 

21) 

L3*63  4A.'«0  EOGc 
3AnJ  SwGE 
0E5!»ED  '*LJ£ 
•cIGht;:**; 


a**o  i 
3,3140900*00 
c,  10030**90 
3,3300000*0 

i  t?0030ci«0a 


9ano  ? 
a,m5625a0 
a.^aofctfoaao 
-3.333333333 
1.200030993 


9AN0  3 
a.49t562«30 
9.320090000 
9,300009090 
1,300900909 


E*TPE«4L  *953UE<CIE3 


0,3000040 

0, 0390*25 

3,3671375 

0,9090625 

2,1  994009 

3,1  109375 

3, 1312500 

0,1562509 

0,1412500 

9,2976125 

9,2324125 

3,2593750 

3,2059375 

0,3125004 

9,3375999 

3,3640625 

3 , 3a  75330 

a.40i562? 

0,4140625 

3,4464759 

a,*aa9ta3 
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i  ‘e,Oc  :r:  ^ 

l;4-* :ro:,'*L  *11.75*  oe$rs* 
^  VLr.C*  !  r-« 

«am)**S5  ','tur£3 

*T Lr6*  LENGTH  a  75 


12-6-  *A-\0  633t 
^ot.^  SA'IO  cCGc 
3*3;5g:  #4^.5 

;£y f i t ; o  s 


»*  f 

n « 

7,150*5**46-01 

* 

"( 

76) 

~r 

25  * 

•Ja!  1  4.7*1  2°€— 10 

a 

*1  f 

7  O 

"  ( 

i)  « 

4,242443516-01 

a 

4( 

73) 

•*  ( 

4]  ■ 

-«l.9202130-6-J2 

a 

"C 

7?) 

**  ( 

55  « 

-.I2«*7a2l£-0l 

a 

-r 

71) 

a)  * 

-0,740445796.02 

a 

7( 

79) 

M 

7)  * 

-0,1911/4476-01 

a 

4( 

*9) 

M  c 

*)• 

0,195314876—11 

a 

*< 

*8) 

-< 

9!  a 

-3,2”204925€-n 

a 

*7) 

-c 

151 » 

3.2-72**576-01 

a 

**) 

*c 

m* 

0,3734i 24?5,o* 

a 

-7 

*5) 

hc 

123* 

-0,223*35156-01 

a 

4f 

*4) 

« ( 

13)« 

0.251 14-7 J6-01 

a 

-C 

*3) 

►*( 

1  -  3  * 

-0,4444197*6.21 

a 

-f 

*2) 

*c 

153* 

0,23770**06.01 

a 

*1) 

*< 

la)  ■ 

0, 12! 0t7*l£-Jl 

a 

*0) 

nc 

I  7 }  * 

-0,2^*304346-01 

a 

-( 

59) 

«( 

193  * 

0.5c62255lf-0t 

a 

58) 

*c 

143  ■ 

-0 ,587283256-0 1 

a 

4  7 

57) 

Hc 

223  * 

2,113474376.21 

a 

56) 

-*< 

an* 

0.333*754*6-02 

a 

-*( 

55) 

*C 

22)  ■ 

-0. 4332575*6.01 

a 

4( 

54) 

-*( 

23  3  * 

0,854322476-01 

a 

-r 

53) 

2fl3* 

-2,442793786-01 

a 

«< 

52) 

*7  ( 

233  * 

0,1853494*6.01 

a 

51) 

"( 

2*3  * 

3.43*173146-02 

a 

501 

«( 

27)* 

-0.922012476-01 

a 

49) 

?5)« 

0.788042246.01 

m 

•9! 

48) 

m  ( 

24)  « 

-0,340934856-01 

a 

*( 

47) 

«i 

301  • 

7,-30112096-01 

a 

-f 

4*) 

313* 

0.72bt7342€—U 

a 

M  f 

45) 

-*c 

32)  ■ 

-0 , 1249*7636*00 

a 

-( 

44) 

333  * 

3.3ca.T4a57?-at 

a 

4  C 

*51 

34)  ■ 

-7,8433299*6.01 

a 

-*  t 

42) 

35)  * 

-0, 3?  5*4*  b*6-2 1 

a 

4  7 

-t: 

-( 

36)  * 

0,32*2*7946*00 

a 

-  c 

*2) 

*c 

37)  • 

-0.404**0276-02 

a 

-»( 

39) 

3  8)  ■ 

0.bO2O53556*0O 

a 

-f 

39) 

9  4NH  1 

3  .033000033 

3.143759733 
l  ,30*330*09 

2,  rusaaaa** 

3. 771 913*69 


9  a  no  * 
3.1*5339430 
3.3733345’ 3 
0.303330020 
1  .2  3^3  33330 
3,  ’33295401 


JE*:*7;g?*  JS  2a  -22.  8  75833879  .9 .  0* 1  8  3*353  - 


9  ANO  3 
3. 416797443 
2.930300333 
1  ,  ,130303000 

3,239532750 

0. J85522*85 
2i .35843-3-* 


£»r*c-a^  »4r^ugNci£3 

3,2131574  0,;>37ljaa  a.?4tli84  0.2553°97  3,3*40789 
2.763<;592  3,2472395  3. 1110197  3,1253030  3.137335? 


3.1-37547 

2  ,  *  4*5904 

J. 26**223 
3r  319  3? jo 
2 . 5  b  7  b  4  3  j 
3, 4597232 


l,ta?3344 
0 . 2  J  9  1  7  2  * 

,-i  5h49<73 

?!32«*t3l 

2,37032-5 

3.-  72  7  ^5 


■}.  1*9*299 
7,2234331 

^2^21157 

7.3-! 3**3 
2,42*7317 
?,43587b4 


3. 17*0322 
'*,2314513 
'1.2*3*294 
J, 352"! 7: 
<?.  4324224 
J  .  3  000003 


2,  1  35943b 
7.2-439*9 

7.3354*a4 

7,3*1 3*31 
*5,3455903 
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TABLE  9 


p;m:*e  i^ulsc  response 

u:nfar  ®wAsg  orcx’Ac  pjlte9  jesio* 

°emse  £xc«*n<;?  algapithh 

RanoPaSS  FtL^C® 


LCfcCR  5A.40  EDGE 
UMfR  9anO  EDGE 
3C3IR53  v* tUE 
mcISHTING 


p: 

L 7C»  LENGTH  a  7? 

•• 

••*  JHRUL 56  PC3P0N5C 

•  • 

•  *« 

fit 

;)*  9,iia779P4C^ai 

»( 

75) 

*c 

ai»  a.aaT«M<»§e-aa 

7  “  ] 

"( 

3)4  0.21593887E-*: 

*•< 

73) 

4)4  -0.37731 763C. 32 

Hf 

72) 

51  *  1.a»in»''8C.« 

71* 

4)*  -9. 1 t0l?09:E-01 

*< 

70) 

HC 

7).  •a,427223]i(.t) 

69) 

ai»  a.tat57i3je-aa 

Hf 

68) 

*C 

9)*  -0. 4*3134535-92 

*( 

67) 

10)*  3.U310730E-31 

66) 

»c 

it)«  fl.a«ia«a*ae.aa 

45) 

*c 

12)"  -0.923ai427E.92 

64) 

HC 

13)*  -0,232840 9fe £-92 

*( 

63) 

*( 

143*  -0, 15925940C-01 

62) 

MX 

15)*  9,2222?040€-*2 

N( 

61) 

H{ 

16)*  3,11 5949ia€ -9 1 

*f 

69) 

H( 

1 71  *  a.4«a77.„E.a3 

59) 

*c 

10)*  3,l2372fl94E-*t 

N( 

56) 

"C 

H)*  -0, 16945603E-01 

N  X 

57) 

20)"  -0. l680tai3E*0l 

*C 

56) 

H  X 

21)*  0, l 0l9t  93JE-02 

*( 

55) 

H{ 

22)*  -0.90078090E-02 

Mf 

54) 

23) ■  0,33937333£-0l 

53) 

*X 

24)*  3.99254601C-02 

*C 

52) 

*( 

25)-  -0. I343aa92£.0t 

*( 

51) 

*C 

26)*  -9, 224369025-92 

Hf 

50) 

27) »  -9.492430575-01 

49) 

MC 

285*  0,  U6378735.01 

48) 

29) »  9,26744b96£.01 

47) 

30)*  9, 134693135-*! 

*( 

46) 

H( 

31)  ■  3.59*242545-01 

*C 

45) 

32)*  -0,700722935-91 

H  f 

44) 

*C 

33)*  -0 , 455609 T 25 -0 1 

*< 

43) 

341*  -0.29541502E.01 

42) 

35)*  -*,647320345.01 

«•( 

4  1) 

9( 

36)*  * , 30« 1 8225E .10 

40) 

*< 

37)«  *, 5697!0l6E-*l 

*  c 

39) 

*C 

36) «  0,521 622 155.09 

"C 

38) 

9AN0  1  9aNQ 

2 

5**0  3 

a. 

000000P30  a, 144582099 

9.159822368 

3* 

1437597J0  3.159000002 

0,1 653*0000 

«• 

000000300  -6t ,022540939 

*1 1 ,666666667 

a* 

045888*90  0,004588009 

0,904589809 

L0«ER  3ANO  £0G£ 
*|P*E0  3 ANQ  E3GE 

ossmo  value 
•CISMTIng 


3ANO  5 

0,37:146073 

9.3d3ep«00« 

6,1326742*5 

3,333527390 


94*0  6 

0.335822 363 
3,43211 3399 
47,897*76037 

3.303853275 


8**0  7 

3,4333422^3 
3,4;6707443 
4.131317733 

9,038532*5* 


IXT0E"*L  M!2 
3*3115132 
3,3333592 
3,1437397 
3, 1973724 
3,2379274 
3,31*4953 
3.37*3243 
*•*579:39 


U6^c:E3 

a, *24315* 

9,3973344 
3,1433224 
3,2345855 
3,2732632 
3,331*43^ 
•7.4029423 
3,^7*896  1 


3.34U  18* 
■*.111*197 
9.1941118 
7,22042:: 
3,2825987 
3.34!  4339 
■’,4:6  797  a 
C, 45587*4 


3,3553987 
3 ,  : 253099 
3,i?65t22 
3,2332546 
9 ,2949342 

3.3541447 

8.4291330 

0,3000040 


3,3490734 

3,1373355 

3,13639:4 

3,245595: 

3,3372697 

3.3443355 
3, 4*39354 


31*0  a 
3,145822348 
3,370324514 
3. 3093*3*09 
l  .300020*03 

•  A*o  3 
3.41*6:9809 
3 .530000300 
3, J02000300 
3. J33532,53 


TABLE  10 


•:*:*£  ;^®ul5e  *E3*:dn3e  c ft*) 

1.1*64*  **43E  DIGITAL  FILTE»  7ESI5H 
^“EZ  exchange  ALG09I7** 

3AN0HS3  FtuTfff 

FIl7£*  LfNG*H  a  75 


•**••  IHPUC3E  *53*0*36 


H( 

Of 

0,1071 44606*9 1 

■ 

HC 

75) 

*( 

25" 

0.32081976E-02 

« 

M  < 

74) 

HC 

33* 

0,l*4l3l80e-9t 

• 

HC 

73) 

"C 

4  )  • 

•0.1 623 t 1 926-92 

a 

HC 

72) 

"c 

53  * 

3.77142632E-02 

• 

"C 

70 

*< 

4)* 

-2, 10*831636-01 

9 

70) 

«( 

7)* 

-9,66070236^.02 

■ 

H( 

6*) 

Hc 

8  3  * 

-9.2*34*4526-02 

a 

H  C 

68) 

"C 

*)■ 

9.17539442E.93 

a 

H  t 

67) 

Hc 

133* 

2 ,060431618*02 

a 

HC 

66) 

*< 

to* 

3,362834448-02 

a 

*C 

65) 

123* 

-9,5279024  0-02 

a 

H( 

64) 

H< 

133* 

•8,71431120-02 

a 

H( 

63) 

HC 

143* 

•9,482434268-22 

a 

HC 

62) 

153* 

-9,1 4444323C *22 

a 

HC 

60 

*c 

163* 

0,12065*66-01 

a 

HC 

60) 

*c 

03* 

1 , 44 1  7 1 104C-02 

a 

HC 

5*) 

*c 

193* 

0, 36*7398  0*22 

• 

HC 

50) 

"C 

1*3* 

-0,349731346-02 

a 

H  ( 

57) 

hc 

233  * 

-0,21640 2* 86-91 

a 

HC 

56) 

HC 

20* 

9,111*47846-02 

a 

H  t 

55) 

hc 

223  * 

-9,36438930-02 

a 

HC 

54) 

H( 

233  * 

8.219*608*6-01 

a 

HC 

53) 

H< 

24)  » 

9,185259026-01 

a 

HC 

52) 

*c 

253  * 

-0,18523**36-01 

• 

HC 

50 

HC 

263  9 

-0,274280096-02 

a 

HC 

50) 

*c 

273  4 

-0,41 *04**06-0 l 

a 

HC 

44) 

hc 

28)  • 

9,532330606-93 

a 

HC 

*8) 

H( 

2*)* 

0,37 9607  4  0-81 

a 

HC 

47) 

hc 

303  • 

9, 5  3053  t  246-22 

a 

*  c 

46) 

*< 

30* 

9,5732*7296-01 

a 

HC 

43) 

Hc 

323  « 

-9,614663466-01 

a 

H{ 

44) 

hc 

333  * 

-0.577222*76-91 

a 

HC 

43) 

*c 

34)a 

-0, *87447166-02 

a 

H  ( 

42) 

Mc 

35)  * 

-0,6*4317236-91 

a 

H( 

40 

*C 

363  * 

3,390739176-00 

a 

HC 

40) 

pi  r 

37)  » 

0.67543830-01 

a 

H  C 

3*) 

*c 

38)  » 

9 , 5  3*2 1 ; 7  66-90 

a 

H  C 

38) 

8AN0  1  3  ano  2  3 ANO  3  3  anQ  a 

3,3380*008*  0,1443820*8  3,16606311  3,37::a%478 

0,143739730  0, 165333*4 3  *,37«32«3i9  a,4ib7*744a 

1,209909939  100.339099090  9,339904090  1 99  *  309030090 

9,344899999  9,949490009  1,390909990  4, <333033490 

3AN0  3  84*0 

U3*£*  3a.no  COGE  3,Al7bt«?i9a 

'jF*E*  9iN0  50GE  3,400090090 

0631*60  ''AiUC  1.390909090 

-CISHTING  9.338832740 


£XT9Ehal  F«e;uE*C:E3 

3,012334$  3,0271382  9, 1*11134  3,3540487  3,36*0734 

3  ,  38  335*2  3,2*733*4  0,  1  1  19147  9,1243000  3,1373355 

0,1437547  9 , 1 66  1  5a  3  9.  16*4459  3,1  768471  3,1947154 

3,1*74063  3.23a*t«5  3.2212543  2,2335*J5  3,244926a 

3 , 25826 1  #  9,2795*71  3.252*326  3,2469405  3,3984260 

3,3257616  9,33227-»7  3,3446123  3,3553919  0.3651645 

3,3733243  3 , 3*55 1 7*  3,4139877  3,4283136  3,4a3ti32 

3,4579433  a,47la*6l  3,43587*-  3,5293099 


L3MC0  8 ANQ  EDGE 
;jF*E9  8AN0  £QG£ 
3631*50  VALUE 
*EJ5ht:*G 


•M 1 1  II  |  II I mU  I 
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TABLE  11 

P  *7*  u*  •:*  u'1  v’»  o  w  »i*  u.  —  >'>j  —  »*•>  •**■>  <*'j  -  V  »7-  • .  ^  ^  t  t  1 

*>  «*s,  '  *  “  -  ■-  *  ^  *  >,  .■*»  •  -  -  ■* 

*  *p  tt  co  —  3  ^  -::  *.  c  p  -r  »***  or;  o  p  **■  f"»  **j  —  3  *t.  r**  k.  \y  p  -r  <**  eg  »-  rz  o  ■> 

,S.  N  N»  ^  <3  <0  **,*  '•L  O  ■  0  J  il‘  'v  \*J  b*j  *p  u  >  «">  P  •*  *  ‘  'P  |  '.  V  T  ~  T  V  *T  T  T  >T  T  C*>  C  S 

» 

WW  y  uWS(  *W  WW^  w  ^  v  ‘<«  V-^  ^  ^  *•  W  %»  VVwvV^/VV^  J 

'•  ■<  ■  >  i  ”  '*  :•  ■'  •'  •*  "  *  ..  i:  ;  ;i  u  'i  ii  .i  **  :i  ?i  •  ,  •, 

i  •**  ~  *\:  ~  'V ^  'i*  ^  «>.<  *i'  z  t  r  t  C-J ^  -1  r  iJ  *  "  -  r  r  r  r  ^  f  x 

_  r  ~  "  r  ~  "'"T"  .  m""  i'mTTT*’"''"  "’""’IT""  :  7 

—  *«i/  tj  ->■*  —i  *.J  ’JJ  'jj  '-J  ijiJ  Uj  *  J  l«*rf  !  '»  *J|  yj  ■. J  ^  'ij  y.  Ck  *jJ  -J  '- J  '-i_>  LlJ  '*i  Ld  '-U  '.u 

—  u-  •*:  u  -  vj  w  **>  o  v  \o  ::  O  —  P  ■::  -  ~  .:  *  :  -  t  -  -  —  •; ;.  •*  CO  p-  -»j 

—  X  u*  •  0  •>  za  r»  ec  —  <’  r*  33-2  —  :  —  *4  r  ,  ;;■  •  •;•  ^  •.*  *•■*  *■•*.  —  •  ■;*  ,*•••  —  -r  7*.  —  *i*  r-.  •  ^  ;o 

'.  j  ^  ■*•*  *r  .“■  v*  (*  <  h.  a  r  o  -•  <■■.*  sd  ^  *r  *7-  •***•  •  j  •'”  *r  *7-  ■*•  •*•*  »t.  .ti,  m  r .  u-.  —  .  %  ,*.j  .jst 

^  j.  ^  •-  •:■*•  t  -■  -* ■-■  r  -  •-• :  -  -r  * >:•  -j  r-  ■/>  u*  ■  •;•  .7*.  :•;•  '•  -  *r*  /•  7  ■«  =*g  •-  ^  t  —  *,T-  r  -  >7.  ■ 

_>  <*g  'L'  »T'  *r  •  i  -  *"  -l*  *0  05  »**•  v!  -r  f  -  '.‘g  •*•*  -.*■  *  -  »'■.' :  .  •*. .  •-  j  *x*  >: j  '*-j  ''*j  *x  -r  *r  *r  4'  * 

—  —  *■“  -  “■  •*  r -  '7i  “  •*■  ^  *7  CO  U"  —  u*  •-  u“  —  — •  — •  *^4  “•  ■  '*•  -i?  ^  u“  -''j  -•  -jj  r.  *r  '*•«;  •*?  a:  4r 

31  H  Sill  IS  I  i  ill  Ills 

*  •  n  .‘  II  I!  il  )l  'I  Ij  M  II  M  (|  l|  H  tl  II  It  II  It  H  It  i»  If  It  H  II  II  II  II  II  (|  -I  ,|  I  'I  ;|  -I 

_  #  |J  '  ■'•  T  u”»  CO  *7‘  *>J  *>7  T  U7  *C»  f-  0j  >T>  ‘3  — ■  |"»J  »*•'  *T  u*  ‘V  ^  CO  ’7'  -*-•  »%J  '■•*'  *T  :;“|  P--  ■!•*> 

X1  ♦  *“  — *  •*  •*+  —  ^  ^  •  -j,  Cg  '*'j  ''»i  **'j  >  u  *u  C'j  '“'J  *****  CO  *  0  i’O  >*•*'  **<  i'O 


w  a*  <7-  uigi*  U*  oj  —  co  •*•*•  <y  r_-  -r.  r-  ‘7*  ~r  I  «  «r  7  o  -  CO  X  *j  ^ o  *t\  <*r  u. 

X  ,vj  °  uo  o c*j *\i oj  —  *3 c* o  ^07- co  o  1*0  -r »>o o  .v  co  —  <M o «w  *r * 0  *‘0 u» it.  ^  ^  x  x;  u. 
•  -3  0  O  03  *3  ’*0  00  o  03  O  O  «.*0  *3  CO  O  O  O  O  CO  >23  ‘3  03  3  O  CO  *33  CO  3  •=•  *3  O  *T»  .jj  .;o  CO  *T  *3  Tv 


JJ  O  -3  o  0  3  O 
,  T  3  <JS  'w  O  O 
~  '"Z>  3  3  w  u"»  ‘4"0 
ZZ  <Xf  CO  7777 

3  IX  iTk  .7S  O  O  .jv  <J> 

X  U . 

r  ct  ii 


ilO  O  3 

M  :0  O 
7 -  03 

•  v  O'  ti3 


S  0  3  3  0  *3  3 

—  F  ^  ?  r  ’='  'X 

X 

if  o  *r  *»r  —  —  3- 

•  >^i  -  -  r  t 

*  cc  . 


o  U3  3 
3  C»i  3 
—  3  O 


'P  P  3 
A4  C,  - 
77  /  '7 


•>7  ■*.  X 


U  3  3  '3  3 


>_•  V - 7  0*0 

v.  n  r  -  oj  co 


—  *  —  _  — 


_  W  '  N.  '  '  —  — ' 

s0  u  •  -  - 
-J  —  —  7  4*i 


:» era  tie 
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TABLE  12 


*:v[T£  :^i*L'l5€  CF!*) 

LI^Sa*  P«4S€  nGtTAt  FTLtEb  a€3t<3* 
3tM^2  £«C*a»jG£  4LGQRIT*h 

OA^cPA^a  FTcrc* 
r  I U r £Q  1-c^UTm  a  43 


•  •••»»  IM»UL5c  9CSP0^3€ 

rt  C 

1)0 

-0. 477239538-82 

• 

Hf 

*C 

2)» 

-9,467 142238*02 

• 

H< 

3)0 

-J  .16323213C-41 

« 

*( 

K  r 

4)  a 

9.11555324C.01 

• 

*( 

*( 

5)» 

0, 4549P205C.02 

• 

*C 

*r 

61* 

-0,244668855-01 

0 

nr 

"( 

7)  * 

-0,63l3«66ie-0t 

• 

*C 

8)  • 

-^,2?l53602€-0i 

0 

h  c 

*■*< 

*)• 

0.84450876E-01 

* 

H( 

13)  * 

3.1 '07671  89C *09 

a 

H( 

*  c 

U)* 

0,3341 7526C-0 1 

a 

*r 

12)* 

•0.686 33 022C-01 

a 

Hf 

*  ( 

13)0 

•0,213623745*06 

a 

*r 

«c 

14)  a 

•0,117733405*30 

a 

*( 

"( 

15)4 

3,18309227C*0O 

a 

*c 

16)« 

3.277l3l46£*00 

a 

17)0 

2,1  342iaa?e*<?(j 

a 

"( 

l  8)  ■ 

-0,634525705-01 

a 

*( 

*c 

!*)• 

-O.37987334t*0a 

a 

*c 

29)0 

-0,432 33 667? *03 

a 

*< 

MC 

2  O  * 

3.209067928*00 

a 

"C 

*c 

22 ^  0 

3,6994ib33E*00 

a 

*( 

‘.Q<«E9  8A*n  £uGc 
34^0  €OGc 
Ocil^CU  vAi,'jc 
€  1 f  I  nG 
0EV!ATI3:< 

o€*urra.9  :*  33 


■JAM  0  » 
?#^j30av*^a* 
0.-*r28933.*9 

1  ,aa«*9**f« 

3. *3377.,  <**3 
-82 ,231902.144 


§4*0  2 
3.227548400 
3.2*9301240 
1  .  .)73030tf09 
3,iHMU7Q 
3  .a0a*73(ji|* 

•66.444901 164 


433 

423 

40 

40) 

34) 

38) 

37) 

36) 

34) 

34) 

33) 

32) 

30 

39) 

24) 

28) 

27) 

26) 

25) 

24) 

23) 

223 

3an0  3 
9,344332890 
2,523090909 
0,300300900 
3,2034^6060 
3,322578738 
•32,323948928 


£xr3£«4L  r9ga'je^CiS3 

2.33500-30  3,3170453 

3,j64®.]23  3,*7280i3 

9,251*462  021444:0 

3.344532§  9.353**55 

0,4468055  J, 4737943 


9.33207*3 
0, 22^5464 
9,2^2644 
8,37200*0 
0,53 00530 


3,0463730 
3,2333994 
0, 237?073 
),299oa  * 


0  »  75 °6 39 l 
*,2384 J21 
0,2^93912 
1,4212373 
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TABLE  13 


wire  ippulJe  response  <ptri 

LX**E*R  PHASE  OX6XTal  PILTER  3«XfN 
RE*EZ  EXCMAN6E  AL5Q»ITi*M 

34*0*433  PILTER 

PIUTE*  LE*6TW  •  *3 


impulse  response 


*c 

U« 

4,U3723«7e-42 

« 

91) 

H< 

2J» 

.3.7*43*3479.42 

• 

HC 

*2) 

H( 

3)  • 

•,47*72t*4e-4a 

• 

H( 

«t) 

H( 

*)• 

•4,23344*131 .33 

* 

H( 

*0) 

HC 

9)* 

3,132939*91.93 

• 

H( 

39) 

«< 

9.*2492t32C.92 

* 

HC 

31) 

»< 

7)* 

•#.132979472-21 

a 

"C 

37) 

H( 

6)» 

••.13979333C.31 

a 

H( 

39) 

H< 

9>* 

9,1*33*4472.41 

a 

*C 

33) 

»< 

10)* 

3,2743*a73C.32 

a 

HC 

3«) 

tl)* 

3,29*4l24«e.92 

a 

H( 

33) 

H< 

12)  « 

3,2*7931412.91 

a 

*< 

ia) 

H( 

131» 

•9.3*49713*2.91 

a 

H{ 

31) 

M< 

!•)• 

.9,4*3139472.41 

a 

H( 

30) 

H( 

133* 

9,2913*7722.91 

a 

H( 

a#) 

H< 

t*)« 

3,437492422.32 

a 

H  ( 

ai) 

171* 

9.1323*4922.91 

a 

HC 

37) 

10)* 

9.1 172(3*92.43 

a 

H( 

a#) 

H( 

19)* 

•3.344973342.31 

a 

HC 

23) 

HC 

20)  * 

.9. 274243242.93 

a 

HC 

2*1 

2D" 

3,224374242.41 

a 

HC 

23) 

*C 

22)  * 

4.3479t9**2.<*9 

a 

HC 

22) 

BAWO  i 

84*0  2 

SA0O  3 

*A*0  a 

LOWER  S4N0  2362 

6,3000*0000 

3.374223799 

a. 11*400*59 

0,122*20455 

UP*2R  94*0  2062 

0. 27360330# 

3,1949999*9 

0,121300*00 

2,170391730 

023X420  VALUE 

3. *74*972*3 

4.100*4**07 

U, 5239*70*2 

•EI8MTIN6 

3. 3 19993999 

0,001000000 

% 

0,0001*0320 

94*0  9 

34*0  * 

0ANQ  7 

8AN0  S 

L0.2R  94*0  2062 

4.1*4772294 

3.1427722P9 

0,2200*4*09 

3,291721**0 

UPPER  *4*0  2062 

3.141391799 

3.22794*439 

0,2903012*0 

3,329000000 

023XRE0  VALUE 

4, 1*****4*7 

9.*7»*973*3 

0,000000000 

•3,346330*78 

wEI6*TI*6 

3.991*93299 

9.91*932949 

0 , 1*0320000 

0,203**0600 

94*0  4 

84*0  13 

0ANO 

LOWER  94*0  E06C 

3,313*29499 

8.343433243 

UPPER  *4*0  2062 

3. 3*4932*33 

2.339999999 

023XRE0  VALUE 

•3*,33**54*?9 

3, 393999993 

•2X6*71*6 

9,4943**9*4 

9. 3334*4893 

e*tre*al  prequenciss 

•,3333333  3,3134333  3,933331*  3.9333773  3.39973t4 
3,373*333  3,1399393  3.1337722  3.23*3972  3,2349347 
9,2«**933  3,2*44*43  3.27**34*  3,2««*2«3  4. 3343399 
9,3499933  3,399*4*4  9.3772333  3.3344399  4.4394337 
9,444*434  9,479214*  9,3933399 
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TABLE  14 


unite  inruloe  resronse  (rim 
LINEAR  PHA8E  OIGXTAL  filter  oesxon 

»[H12  exchange  ALQORITHH 
BAN0FA99  FILTER 
FILTER  LENGTH  •  *3 


8**0  3  8**0 
LONE*  SAND  BD«E  8. 3*3333333 

UFFER  8*N0  EDGE  9.988989089 

DESIREO  VALUE  8.08*889808 

NEIGMTING  8.983*88888 


EXT*E*»L  EREQUENctEE 

8.8809908  8.03)3088  8,9*38138  8.38I0T93  9,8739033 
8.1133781  8.1*39381  9,1733331  0.3998**3  9.3389893 
8.3*3173*  0.3893193  8.37*4339  0.383797t  9,3903013 
8.3*39333  9.333938*  0.3737839  8.3999383  0, *39*997 
9.*«98*9*  8, *733148  8,3008000 


i  KHuni  45  tmtiipAvs  nt.ie* 
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TABLE  15 

o  CU  »7%  \j>  A  V:  Z  U.  'll'  OW  w  -•  Ui  it*  CD  00  CO  U  «  Up 

£ 

^  O  O  — •  U"'  'VI  N  ••  ••  'ii  t3  O  to  li*)  u”  CO  |V,  00  |£i 

J*  CO  0!'  ■=■  CO  o  O  CO  'X*  ®  ®  ■-  CO  C*  »*'J  *7.  Ui  ®  f  - 


O  C*J  (J-  IS  ••  UVI  L  ’J)  ' J  00  u  ••  T  U  U"  CD0>  CD  Ui  O  L. 

w 

3:" 

fcOOPO^’S’OCOODOOOO^^OO  ©»\|  cr.  UJON 


Copy  available  to  DTIC  does  not 
P*mut  fully  legible  reproduction 
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TABLE  16 


c  *C3*0*Je  C*t*> 

LI'£A*  .jISTTac  *TLTr«  0€3I*N 

>t*»Z  F^ChANIJC  \LC,Q*X*4* 

m»i>0*22  rTur2* 
rtUT !•  LENGTH  «  57 

•  •*••  I*®Ul5€  *E$*nN3E  «••»• 

«<  n*  -0.4403tS2t«-*3  •  Hf  57) 

-(  2!»  -n#*3«33*a4e.42  ■  ■*(  34) 

n(  3) «  -0,447243410-02  •  *(  ss) 

*<  4)*  -«,t520462*«-0i  ■  «•(  3«) 

■*(  3)  ■  -4,1  *0733431-01  ■  *(  S3) 

-»(  *)•  -4,t732«942«-»l  •  H(  32) 

H(  7) »  -0,J3l2799«0-0t  •  H(  51) 

«c  a)«  -*.t34j72«3e-*2  *  *t  30) 

Hf  4)«  a,635732*9C-03  »  *t  44) 

H(  IJ)*  9,179041)166-01  »  Hf  44) 

H(  U)«  et4774«dl5C-4l  ■  H(  47) 

*C  12)*  9.59440959C-01  •  H(  44) 

H(  13)4  *,l33a3466e+03  •  *<  4$) 

«(  14)*  0'477*44Me-3l  4  H(  44) 

HC  15)4  e\7334<*3«7E-Jl  4  H(  43) 

hc  I*)*  tt,a4aa3s^4€-ai  *  hc  «2) 

m(  17)4  3,31*127016-01  4  h(  4i) 

«(  l •) ®  -*,StS43437«-0t  *  H(  45) 

H(  19)*  -*,14«M74e#«a  •  Hf  34) 

H(  25) •  -*,l22279«9E»tfe  4  *f  IS) 

*c  21)4  «4,2226430SC*00  4  n(  37) 

N(  22)4  -0,993026236-01  4  H(  34) 

H(  23)4  -a,il34437iE#55  •  H(  35) 

H(  24)4  -«,l79?ft*44e«4e  *  H(  34) 

*(  25)4  a,25444645C*5«  *  H(  33) 

H (  26)4  *,19642734C«09  *  *C  32) 

Hf  27)4  9,22394636E«00  4  H(  31) 

H(  261 •  »,14S**aTaC*«*  *  h(  32) 

H(  2«)4  2,7as«4l73€*29  *  n<  09) 

<4A»o  i  8a*0  2  6  ano  3  Sanq  a 

4ano  ec6€  0,11*169330  2,244313732  2.374735312 

’****«  rU*»o  EJGC  2.72422222^  0.2*3034473  *.311367133  2,349244562 

G£37s€u  VAU.'c  1.7075^7^7  2.330909000  1,00*9*9*90  2,200404495 

.erSnrihtt  2.143265262  9. 17459*272  1,249259000  3.141602569 

:evlAUON  0, *0*347446  0,3996l42U  9,30*156136  0,600394644 

u€^r*rl3N  jn  09  -65.372632U4  -44, 1632301 1 1  -74,320631374  -46,514632221 

6  amO  3  SA*0 

L'JwCR  3  A  NO  60G6  3, 460991 790 

00*59  3A 40  EUGC  3,509090909 

C53I9£D  Vii.u€  1,90.1090094 

u6IG«7InG  0.21 3 l 96449 

wCVlAfXON  3.002597142 
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9,1191604 
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9.1306423 
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9,1534470 

2,1623664 

2,1942401 

3,2336347 

a.aMjtjr 

3,2714464 

0.2799123 

a.aaartaJ 

9,2994641 

a.3a«ts6« 

3,3115673 

0.3747355 

J.JTMUf 

0,3633362 

7.44299! 7 

t?, 4602045 

9 «  *077742 

a(sa««.a. 
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"t 
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*( 
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»( 
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"< 
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"f 
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»< 

34) 

N( 
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"f 
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*< 
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*< 
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2*)*  •*. 119*193*6.01 

*< 
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H( 
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9.31*3400*9 
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003334*41** 
0.341 34037* 
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UMER  0*w0  6066 
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*7,144437143 
3,417494939 

94*0  * 
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(*AM0  7 

0. 1 t  t  4*3436 
*,203034*70 
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§4*0  4 
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3.310*42**3 
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3.21300030* 
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<4*0  CQCC 
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S*NO  13 
*,3199*42** 
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*, 0*129*6*9 
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9.37047790* 
9.37*735513 
.*.47595*9*4 
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5AN0  13 
2.377*1394* 
1.3**244300 
*.*00200000 
0,1*1*02980 

%A"0  t% 
2.34*37794* 
3,49*3*0303 
0,33*723634 
3,316160254 

W3*C*  IlNO  cost 

UMC*  3**0  COG C 

oesifffo  V4cue 
«<X5*rX«G 

04NO  17 
•,*9957791* 

0, **77412*0 
22.1221*413* 
*,321314**4 

9AN0  10 
l.**«0*092* 
*.**2441X4* 
9.344723950 
2.21314*440 

4»NO  19 
2.4949*432* 
2, {0*409400 
2,200004000 
9,21314*440 

3an0 

6*T*6“*L  *966U6HC169 

9. 24,2909 

0. 9,020*2 
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*, *4 1390* 

». 774,3776 

9.0999310 

»,U773l{ 

*.l3t 7»*l 

*.14*42*3 

3,1*44*01 

a.l790T<k] 

0,i9«a*0i 

2.205*597 

*. 21*090* 

3.371446* 

I,*790t20 
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9,2994**1 
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HC 
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a 

HC 
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a 
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H( 
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a 

H( 
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»( 
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4 
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M( 
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*< 
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M( 
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HC 
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a 
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HC 
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9 
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HC 
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9 

HC 
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H< 
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*C 
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a 

H< 
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a 

HC 
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we 
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a 
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«c 
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a 

HC 

«4) 

hc 
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a 

HC 

34) 
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a 

HC 
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H( 
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a 

HC 
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a 

HC 

34) 

Hc 

23)7  (.Sl*7«77*e.8l 

a 

HC 
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Hc 
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a 

H( 
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H( 
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a 

HC 

33) 
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a 
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34) 

Mc 
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H( 
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L3»(ft  8*N0  (066 
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9.11127973* 
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4,494492444 

9.113139399 
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0(919(0  V*UU( 

l ,194419444 
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199.399999999 
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9,173993299 

9,317999899 
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4.417444444 

4ANO  9 

SANO  A 

9AN0  7 

BANO  4 

LOrtCff  SAMO  C04C 

9,27934)314 

4.312444734 

4,377413944 

9,3793771** 

UMCft  (4*0  S066 

4,311347199 

9.37*733919 

4.4*244)744 

0(317(0  V*LU( 

198.339999888 

194.449444444 
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8**0  7 

8**0 
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9, 4**997329 

U77(«  8**0  (062 

9,999988979 

0(317(0  7*UU( 
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N(X6hTXN6 
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(me*HC  MCSUtNCin 


9,3*92922 

9,93*9*43 

3,999*1*9 

9,9779792 

9, It  12*9* 

9,1199979 

9.l3t7«91 

9,1*77*37 

9,1*277*1 

9,17997*3 

9,17297*7 
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9.2939977 
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9,39*19** 
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3,34173** 

9,37787*7 

9,3*7377* 

9,43*3*99 

9,4297379 

3,4*7*973 

3,493**97 

9.39*9999 
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"t 
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h( 
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*c 
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9.34.774479*49 

a 

*2) 

133* 
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*c 
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fit 
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fit 
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*c 
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fit 

52) 
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fit 

513 

i 4)  a 
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*C 
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